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NON-GAUSSIAN  NOISE  MODELS  AND  COHERENT  DETECTION 

OF  RADAR  TARGETS 


1.  INTRODUCTION 


The  detection  of  radar  targets  agaiast  a  background  of  unwanted  clutter  due  to  echoes  from  the 
sea,  land,  or  weather  is  a  problem  of  lundamental  interest  in  the  radar  community.  Attacking  the  prob¬ 
lem  generally  requires  an  understanding  of  the  statistics  of  clutter  echoes.  Insight  into  the  possible 
structure  of  the  optimal  (or  quasi-optimal)  detector  can  be  gained  only  after  the  statistics  are  under¬ 
stood.  Hence,  the  task  of  formulating  optimal  detection  processors  separates  naturally  into  two  funda¬ 
mental  steps:  1)  formulation  of  a  statistical  model  for  the  noise  from  which  the  desired  target  is  to  be 
discriminated,  and  2)  application  of  the  results  of  detection  theory  to  determine  the  structure  of  the 
detector  for  the  type  of  noise  modelled  in  step  1.  (Although  noise  often  refers  to  the  ubiquitous  ther¬ 
mal  noise  in  the  receiving  equipment,  here  it  refers  to  any  unwanted  background  disturbance,  including 
clutter  echoes.) 

For  many  years,  radar  systems  had  relatively  low  resolution  capabilities  such  that  clutter  echoes 
were  thought  to  comprise  a  sum  of  responses  from  a  very  large  number  of  elementary  scatterers.  From 
the  Central  Limit  Theorem  (CLT)  of  probability  theory,  workers  in  the  field  were  led  to  conclude  that 
the  appropriate  statistical  model  for  clutter  was  the  Gaussian  model.  With  this  model  in  hand, 
researchers  were  able  to  attack  and  solve  the  detection  problem.  The  optimal  detection  structure  was 
discovered  to  be  the  now  well-known  matched  filter.  The  importance  of  this  result  for  the  radar  detec¬ 
tion  problem  should  not  be  underestimated.  Tlie  structure  of  the  matched  filter  may  be  viewed  as  that 
of  a  whitening  filter  followed  by  a  correlator  that  maximizes  the  output  signal-to-noise  ratio.  The 
structure  of  the  frequently  used  coherent  detection  processor,  which  is  composed  of  a  moving  target 
indication  (MTI)  filter  followed  by  a  bank  of  Doppler  filters,  is  essentially  a  suboptimal  implementation 
of  a  bank  of  matched  filters.  The  MTI  filter  fulfills  the  role  of  the  whitening  filter,  whereas  the  bank 
of  Doppler  filters  plays  the  role  of  a  bank  of  correlators.  In  a  sense,  this  coherent  detection  processor 
may  be  said  to  have  been  derived  according  to  the  two-step  approach  described  earlier. 

In  the  quest  for  better  performance,  the  resolution  capabilities  of  radar  systems  have  steadily 
improved.  For  detection  performance,  the  belief  originally  was  that  a  higher  resolution  radar  system 
would  intercept  less  clutter  than  a  lower  resolution  system,  thereby  increasing  detection  performance. 
However,  as  resolution  has  increased,  the  statistics  of  the  noise  have  no  longer  been  observed  to  be 
Gaussian,  and  the  detection  performance  of  the  traditional  detection  processor  has  not  necessarily 
improved.  The  radar  system  is  now  plagued  by  target-like  "spikes"  that  give  rise  to  non-Gaussian 
observations.  These  spikes  are  passed  by  the  detection  processor  as  targets  at  a  much  higher  false 
alarm  rate  than  the  system  is  designed  lo  tolerate.  The  mason  for  ilic  poor  performance  can  be  traced 
to  the  fact  that  the  traditional  radar  detection  processor  is  designed  to  operate  against  Gaussian  noise. 
New  detection  processors  are  required  to  reduce  the  effects  of  the  spikes  and  to  improve  detection 
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performance.  According  to  the  two-step  paradigm  described  above,  the  detection  problem  cannot  be 
addressed,  though,  until  the  statistical  nature  of  the  non-Gaussian  statistics  is  understood.  These  prob¬ 
lems,  i.e.,  the  develojxnent  of  a  model  for  the  non-Gaussian  statistics  and  the  subsequent  development 
of  a  detection  processor,  are  the  focus  of  this  report. 

Although  at  first  glance  the  iq)proach  to  solving  these  problems  appears  to  be  straightforward,  the 
solution  to  the  problem  in  step  1,  i.e.,  determining  an  ai^ropriate  statistical  model,  is  quite  difficult 
This  problem  is  complicated  by  the  fact  that  it  is  not  sufficiently  well-defined  to  ensure  a  unique  solu¬ 
tion.  As  stated,  the  goal  is  to  determine  a  model  for  non-Gaussian  statistics,  yet  this  problem  undoubt¬ 
edly  has  an  infinity  of  solutions.  Which  of  these  solutions  is  proper  for  the  radar  problem?  To  answer 
this  question,  one  could  attempt  to  measure  the  statistics  by  fitting  empirical  models  to  data  collected 
with  high-resolution  radar  systems.  This  approach  is  actually  quite  reasonable  for  determining  first- 
order  statistics  of  the  non-Gaussian  clutter,  i.e.,  the  amplitude  probability  density  function  (PDF),  and 
much  woilc  has  been  devoted  to  this  effort.  However,  two  problems  arise.  First,  if  one  group  of 
experimentalists  measures  Weibull  statistics,  and  a  second  group  measures  K  statistics,  which  results 
should  be  accepted  as  correct?  Second,  the  pulse-to-pulse  clutter  returns  are  usually  correlated  with 
one  another  and,  hence,  require  a  multidimensional  model  for  the  proper  description  of  their  statistics. 
Measurement  of  multidimensional  statistics  is  prohibitively  difficult.  Moreover,  even  if  the  measure¬ 
ment  of  multidimensional  statistics  were  feasible,  very  few  closed-form  expressions  for  non-Gaussian 
multidimensional  PDF’s  of  correlated  random  variables  are  available  to  describe  the  measurements. 
For  those  models  that  are  available,  no  useful  criteria  for  deciding  among  them  are  known.  Note  that 
this  latter  problem  also  arises  for  Gaussian  statistics.  Even  if  the  first-order  statistics  are  Gaussian, 
modelling  the  multidimensional  statistics  by  the  multivariate  Gaussian  PDF  must  be  justified  in  any 
particular  problem.  Fortunately,  as  indicated  earlier,  the  choice  of  the  multivariate  Gaussian  PDF  in 
the  radar  problem  aj^ars  to  be  justified  by  die  multidimensional  form  of  the  CLT.  An  analogous 
theorem  is  needed  in  the  iK>n-Gaussian  case. 

The  appropriate  multidimensional  non-Gaussian  model  for  use  in  radar  detection  studies  must 
thus  incorporate  the  following  features:  a)  it  must  account  for  the  measured  first-order  statistics,  b)  it 
must  incorporate  pulse-to-pulse  correlation  between  data  samples,  and  c)  it  must  be  chosen  according 
to  some  criterion  that  clearly  distinguishes  it  from  the  multitude  of  multidimensional  non-Gaussian 
models  satisfying  a)  and  b).  This  problem  forms  the  focus  of  Sections  2  through  S  of  this  repon. 

The  resolution  of  the  problems  listed  above  is  found  in  an  extension  of  the  CLT  to  incorporate 
number  fluctuations,  i.e.,  fluctuations  in  the  instantaneous  number  of  elementary  scatterers  contributing 
to  the  scattered  echo.  This  result  gives  rise  to  a  class  of  first-order  statistics  sufficiently  general  to 
encompass  the  Rayleigh,  K,  and  Weibull  distributions,  as  well  as  mixtures  of  these  distributions. 
These  distributions  are  among  the  most  prevalent  distributions  used  in  empirical  studies.  The  result 
also  allows  for  pulse-to-pulse  correlation  between  the  samples.  Most  importantly,  since  the  result 
arises  from  a  reexamination  of  the  scattering  problem  itself,  it  gives  a  natural  criterion  for  selecting  the 
appropriate  multidimensional  model. 

The  availability  of  a  multidimensional  model  for  the  non-Gaussian  noise  allows  one  to  attack  the 
problem  of  detecting  radar  targets  against  such  noise.  This  problem  is  the  subject  of  Section  6.  The 
goal  of  thi.>  aspect  of  the  problem  is  an  understanding  of  the  structure  of  the  optimal  (and  possibly 
suboptimal)  processor  for  detecting  radar  targets.  Although  the  structure  is  described  by  the  likelihood 
ratio  of  classical  detection  theory,  one  seeks  greater  insight  into  this  structure,  with  the  ultimate  goal  of 
deriving  a  detection  processor  that  is  implementable  yet  close  to  optimal  in  the  same  sense  that  the 
current  processor  based  on  the  MTI  filter  followed  by  the  Doppler  filter  bank  is  an  implementable  but 
suboptimal  version  of  a  bank  of  matched  filters. 
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Section  2  presents  a  general  limit  theorem  for  a  sum  of  a  random  number  of  random  vectors. 
This  theorem  represents  an  extension  of  the  classical  CLT  and  as  such  may  have  application  in  areas 
other  than  the  problem  of  interest  here.  An  essential  feature  of  this  new  result  is  that  it  incorporates 
fluctuations  in  the  number  of  random  vectors  being  summed;  this  feature  gives  rise  to  limiting  statistics 
that  may  be  non-Gaussian  in  nature.  Along  with  this  theorem,  the  notion  of  a  canonical  model  for  the 
probability  mass  function  (PMF)  describing  the  number  fluctuations  is  also  introduced. 

Section  3  then  applies  the  new  theorem  to  the  study  of  scattering  of  electromagnetic  waves.  This 
woric  leads  to  the  idea  of  modelling  the  radar  clutter  amplitude  statistics  by  the  class  of  Rayleigh  mix¬ 
ture  distributions.  Previous  investigators  have  modelled  amplitude  statistics  by  special  cases  of  Ray¬ 
leigh  mixtures,  and  a  brief  review  of  their  work  is  presented.  Because  identifying  a  given  distribution 
as  a  Rayleigh  mixture  distribution  is  a  problem  of  interest,  necessary  and  sufflcient  conditions  for  a 
given  cumulative  distribution  function  (CDF)  to  be  a  Rayleigh  mixture  are  also  presented  in  this  sec¬ 
tion.  These  conditions  are  then  used  to  show  several  examples  of  Rayleigh  mixture  distributions. 
These  examples  include  the  WeibuU  distribution  for  the  range  of  its  parameters  that  are  generally  of 
interest  in  the  radar  problem.  This  demonstration  is  important  because  it  offers  a  physical  interpreta¬ 
tion  for  the  appearance  of  certain  distributions,  such  as  the  Weibull  distribution,  in  radar  clutter  studies. 
These  distributions  appear  to  be  justifiable  otherwise  only  on  empirical  grounds.  It  also  unifies  the 
approach  to  modelling  amplitude  statistics  and  clarifies  the  distinction  between  modelling  clutter  by  the 
K  distribution  or  the  Weibull  distribution;  both  models  arise  from  the  same  underlying  physical  picture 
with  the  difference  between  them  occurring  because  they  have  different  underlying  number  fluctuation 
models. 

The  Rayleigh  mixture  distribution  is  defined  by  a  certain  CDF,  X  >  0,  and  the  determina¬ 
tion  of  this  CDF  both  analytically  and  empirically  is  important  for  the  successful  application  of  Ray¬ 
leigh  mixture  distributions  to  radar  studies.  Section  4  presents  results  on  both  of  these  problems.  A 
particularly  interesting  result  is  that  the  K  distribution,  which  has  been  applied  to  empirical  fits  of  radar 
clutter  data  with  success,  a{^ars  as  the  first  term  in  a  particular  type  of  infinite  series  expansion  of  the 
general  Rayleigh  mixture  distribution.  This  result  would  seem  both  to  explain  the  appearance  of  the  K 
distribution  in  empirical  studies  of  scattering  problems  and  to  offer  an  approach  to  achieving  a  data  fit 
better  than  the  fit  of  the  K  distribution  by  including  more  terms  in  the  series.  Several  examples  illus¬ 
trating  this  latter  idea  are  presented. 

Because  many  radar  scenarios  involve  making  detection  decisions  based  on  multiple  pulse  data 
that  is  correlated  from  pulse  to  pulse,  the  development  of  a  theoi^  of  target  detection  for  these  types  of 
problems  is  necessary.  The  development  of  such  a  theory  requires  models  for  the  multidimensional 
statistics  of  correlated,  non-Gaussian  clutter,  the  formulation  of  such  a  multidimensional  model  is  the 
focus  of  Section  5.  Based  on  a  physical  interpretation  of  the  theorem  presented  in  Section  2,  a  mul¬ 
tidimensional  model  in  the  form  of  a  Gaussian  mixture  is  presented.  This  model  can  incorporate  first- 
order  statistical  information  about  the  amplitude  of  the  clutter  as  well  as  information  about  the  correla¬ 
tion  between  pulses.  It  also  has  the  advantage  of  being  motivated  by  a  physical  picture  and,  hence,  is 
not  an  arbitrary  construction  of  a  multidimensional  model.  However,  the  model  is  based  on  an  ideal¬ 
ized  physical  picture,  and  the  problem  of  incorporating  a  more  realistic  physical  picture  is  briefly  dis¬ 
cussed.  This  particular  problem  is  quite  difficult  and  will  form  the  focus  of  future  research. 

Based  on  the  model  for  the  multidimensional  statistics  of  the  correlated  clutter.  Section  6  explores 
the  problem  of  detecting  a  signal  against  a  background  of  such  clutter.  Previous  work  on  this  problem 
by  the  present  authors  is  briefly  reviewed.  The  problem  of  detecting  a  signal  of  known  Doppler  struc¬ 
ture  but  unknown  complex  amplitude  is  then  attacked  using  the  generalized  likelihood  methodology. 
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The  Structure  of  the  detector  is  shown  to  be  th^  of  a  square  law  detector  whose  output  is  processed 
through  a  nonlinearity.  The  structure  may  also  be  interpreted  as  a  matched  filter  compared  to  a  data- 
d^ndent  threshold.  Section  6  closes  by  presenting  performance  results  for  several  cases  of  interest 

Although  the  report  is  long  and  heavily  oriented  towards  mathematics  and  statistics,  the  results 
are  rxmetheless  usefiil  to  the  average  radar  engineer  who  may  not  have  an  interest  in  the  mathematics 
but  who  is  interested  in  either  the  modelling  of  non-Gaussian  statistics  or  the  detection  of  targets  in 
non-Gaussian  rx>ise.  For  the  worker  interested  in  applying  these  results  but  not  interested  in  the 
theory.  Section  4.2  presents  a  methodology  for  fitting  data  to  a  Rayleigh  mixture  distribution,  whereas 
Section  6.4  provides  both  closed  form  expressions  and  curves  for  the  detection  performance  of  the 
detection  processor  derived  in  this  report  For  the  worker  interested  solely  in  the  problem  of  modelling 
non-Gaussian  statistics  in  scattering  proUems,  Section  2  is  particularly  pertinent  The  results  of  Sec¬ 
tion  2  also  suggest  that  a  possible  physical  mechanism  for  non-Gaussian  statistics  is  a  fluctuation  in  the 
number  of  scatterers  contributing  to  the  scattered  field.  This  observation  in  turn  suggests  that  dynami¬ 
cal  models  that  could  give  rise  to  these  fluctuations  be  studied  further.  Furthermore,  as  discussed  in 
Section  6,  the  availabilty  of  such  dynamical  models  would  engender  further  results  in  the  multidimen¬ 
sional  modelling  problem.  For  the  worker  interested  in  the  structure  of  the  detection  processor  derived 
herein,  both  Sections  6.1  and  6.3  give  simple  interpretations  of  this  detection  structure;  these  interpreta¬ 
tions  could  form  the  starting  points  for  deriving  implementable  structures  that  would  be  advances  over 
the  current  MTl/Doppler  filter  bank  processors. 

2.  A  PHENOMENOLOGICAL  MODEL  FOR  SCATTERING 

In  many  propagation  and  scattering  problems,  the  electric  field  at  a  given  point  may  be 
represented  by  the  phasor  relationship 


£(r ,  t)  =  (1) 

1  =  1 

where  0)  is  the  carrier  frequency  of  the  incident  radiation,  r)  is  a  form  factor  that  determines  the 
angular  distribution  of  the  radiation  from  the  rth  scatterer  and  <t>(r,  r)  is  a  phase  factor  for  the  rth 
scatterer.  Generally,  the  individual  contributions  to  the  sum  are  assumed  to  be  statistically  indepen¬ 
dent.  For  example,  this  kind  of  expression  forms  the  starting  point  for  Chapter  7  of  Beckmann  and 
Spizzichino  (1963)  on  the  probability  distribution  function  of  a  wave  scattered  from  a  rough  surface. 
This  type  of  expression  is  also  implied  in  the  work  of  Dashen  (1979)  and  Flatte  et  al.  (1987)  in  which 
propagation  through  random  media  is  treated  via  the  path  integral  formalism;  Eq.  (1)  would  be  valid  in 
the  so-called  saturation  regime.  If  the  physics  of  the  problem  indicates  that  N  is  "large,"  then  often 
the  CLT  is  invoked  to  conclude  that  E  has  Gaussian  statistics.  However,  in  many  cases,  the  physics 
of  the  problem  also  indicates  that  N  fluctuates,  a  case  to  which  the  CLT  does  not  apply.  For  example, 
in  the  work  cited  above,  Dashen  states  that  in  the  saturation  regime,  N  fluctuates.  As  the  notion  of 
large  N  is  nreaningless  if  N  is  a  random  variable,  a  parameter  that  may  become  large  is  the  mean  of 
iV,  denoted  N.  However,  as  indicated  in  the  work  of  Jakeman  and  Pusey  (1978),  the  statistics  of  £  as 
— >  oo  need  not  be  Gaussian.  They  studied  a  modification  of  Eq.  (1),  namely 

E\r,  t)  =  ^  2  (2) 

yN  i  =  i 
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in  which  iV  has  a  negative  binomial  distribution.  (In  the  standard  development  of  the  CLT.  Eq.  (1)  is 
normalized  by  1/W.  The  normalizations  there  and  in  Eq.  (2)  are  necessary  to  ensure  that  the  mean 
intensity  remains  finite.)  Under  the  assumption  that  the  fl,  ’s  and  <j),  ’s  are  mutually  independent,  the 
Oi's  independent  and  identically  distributed  (i.i.d.).  and  the  ({>,-’8  i.i.d.  and  uniformly  distributed  on 
[0,  27t],  they  showed  that  the  distribution  of  A  =  l£'l  as  ^  — >  «>o,  is  given  by  the  K  distribution,  not 
the  Rayleigh  distribution  predicted  by  the  CLT.  On  the  other  hand,  they  also  performed  this  calcula¬ 
tion  with  N  being  a  Poisson  random  variable  instead  of  a  negative  binomial  random  variable  and 
showed  that  the  asymptotic  statistics  are  Rayleigh.  Thus,  inclusion  of  number  fluctuations  enriches  the 
problem  considerably,  and  a  general  limit  theorem  appropriate  to  this  problem  is  needed.  A  first  step 
towards  such  a  theorem  is  given  below. 

A  conclusion  that  may  be  drawn  from  the  theorem  given  here  is  that  the  distribution  of  amplitude 
statistics  predicted  by  this  type  of  random  walk  model  is  a  member  of  the  general  class  of  Rayleigh 
mixture  distributions.  The  K  distribution  studied  by  Jakeman  and  Pusey  (1978)  is  a  member  of  this 
class.  It  is  also  shown  below  that  the  WeibuU,  gamma,  and  Nakagami-m  distributions  are  also 
members  of  this  class  for  suitably  restricted  domains  of  their  respective  parameters.  These  results  sug¬ 
gest  that  the  general  Rayleigh  mixture  distribution  be  considered  a  model  for  amplitude  statistics  of 
scattering  problems  and  that  the  role  of  number  fluctuations  in  these  problems  be  further  studied. 

2.1  A  Limit  Theorem  for  a  Random  Sum  of  Random  Vectors 

As  is  well-known,  if  Xi,  X2,  are  i.i.d.  random  variables  with  zero  mean  and  finite 

variance  ,  then  the  random  variable 

converges  in  distribution,  as  >  <»,  to  a  Gaussian  random  variable  Y  with  zero  mean  arxl  variarKe 
a^.  If,  however,  N  is  also  considered  to  be  a  random  variable,  the  classical  results  are  no  longer 
applicable.  Theorems  for  random  N  in  which  convergence  to  a  Gaussian  random  variable  occurs 
(hence,  these  theorems  are  called  CLT’s)  are  presented  by  Qiow  and  Teicher  (1978)  and  the  references 
therein.  In  this  report,  the  problem  is  defined  in  a  natural  way,  and  it  is  shown  that  convergence  to  a 
non-Gaussian  random  variable  can  occur. 

Let  be  a  random  variable  defined  on  the  nonnegative  integers,  and  let  (p(n;  N),  N  >  0)  be  a 
collection  of  PMF’s  associated  with  N  such  that  E[N]=N,  where  £[N]  is  the  expected  value  of  N. 
Let  ^  >  1  be  an  integer  and  let  Xj,  X2,  .  .  .  be  a  sequence  of  i.i.d.  random  vectors  defined  in  /?* 
with  zero  mean  and  finite  second-order  moments  given  by  the  matrix  £[X,  Xj]  =  Q  for  all  i  (where 
T  denotes  transpose.)  For  each  N  >  0,  define  a  random  vector  by 
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where  an  empty  sum  (i.e.,  =  0)  is  set  equal  to  zero.  Finally,  for  each  N  >0,  define  a  function  by 

Giz;N)=  id-  j=rpin;  N),  O^z^N.  (5) 

n  =0  ^ 

The  following  theorem  may  now  be  stated: 

THEOREM:  If  G(z)  =  lira  G(z;  N)  conveiges  for  z  e  [0,  <»),  then  converges  in  distribu- 
_  » «« 

tion  as  N  «o  to  a  random  vector  Y  with  characteristic  function  given  by 

Cy(u)  =  Je  2  (6) 

0 

where  F,(r)  is  a  CDF  of  a  positive  random  variable.  Moreover,  Ff(t)  may  be  obtained  by  Laplace 
inversion  of  G  (z ) ,  i.e.. 


G(z)  =  fe-^dF,(t).  (7) 

0 

PROOF:  The  proof  is  presented  in  two  parts.  First,  under  the  conditions  of  the  theorem,  the 
characteristic  function  of  Y  is  shown  to  be  G  (u^Qu  /  2).  G  (z)  is  then  shown  to  have  the  represen¬ 
tation  given  in  Eq.  (7),  which  will  complete  the  proof. 

Conditioned  on  =  n ,  the  characteristic  function  of  Sjj  is  given  by 


C^(u  \  N  =  n)  =  £[exp  \  N  =  n] 


=  (c*(;^)r.  (8) 

where  the  i.i.d.  property  of  the  X,  ’s  has  been  used  and  Cx  (u )  is  the  common  characteristic  function  of 
the  Xi ’s.  From  the  following  asymptotic  expansion  as  lu  l->  0, 

Cxiu)  =  1  -  “  (9) 

(cf  Cramer  1970),  the  unconditional  characteristic  function  of  5^  becomes 
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Cfi(.u)=  p{n;N),  as^->0. 


^.U^Qu  £T  ,u^u.  rr.  U^Qu 

C( — ^ - Noi-^);N)  ifO^ — ^ oi—^)^N. 


(10) 


From  the  definition  of  G(z;  N)  it  follows  that  for  0  ^  Zj  ^  Z2  ^  N,  then  G(z2;  N)  S  CCzj;  N). 
Fix  u .  Since  No  (u^u  /  iV)— >  0  as  N  -*  oo,  then  given  e  >  0,  a  number  K  e  R  exists  such  that 


“Ifi“  .  e  <  “I2“  -  Ato<4‘^)<  +  e 

2  2  N  2 


(11) 


for  all  ^  K,  from  which  it  follows 


Q^u!^  +  e  ;  ^)  <  -No{^^) ;  N)  <  -£  ;N). 

2  L  N  2 


(12) 


Taking  the  limit  as  N— »  «<>  yields 


2  2  N  2 


(13) 


Finally,  since  e  was  arbitrary,  then  by  Eqs,  (10)  and  (13), 


Cy(«)=  \imG{^^^-No{^)-,N) 
2  N 


=  C(i^), 


u^Qu 

2 


< 


(14) 


This  completes  the  first  part  of  the  proof. 

The  second  part  of  the  proof  follows  from  theorem  3  of  Schoenberg  (1938),  which  states  that 
C (u)  =  g(\u\^)  is  the  characteristic  function  of  a  radially  symmetric  random  vector  in  /?*  for  all 
integers  ^  ^  1  if  and  only  if  g(r)  is  a  completely  monotonic  (c.m.)  function  for  r  ^  0.  A  function 
g{t)  is  said  to  be  c.m.  for  t  >  0  if  (-l)^g^^\0  >  0,  /  =  0,  1,  2,  .  .  .  .  The  function  git)  is  furth¬ 
ermore  said  to  be  c.m.  for  r  >  0  if  it  is  c.m.  for  r  >0  and  giO)  =  g  (0^)  <  «».  To  apply  this  result 
here,  let  G  =21  in  Eq.  (14).  This  choice  o£Q  Jias  no  effect  on  the  structure  of  the  function  G(z), 
which  depends  only  on  the  collection  ipin;  N),N  >0).  With  this  choice,  Eq.  (14)  becomes 


C(m)  =  G(Im|2), 


(15) 
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which  shows  that  C(u)  is  a  radially  symmetric  characteristic  function  in  ,  k  >  1.  Hence,  by 
Sdioenbei:g’s  theorem,  G  (z)  is  a  c.m.  function  for  r  ^  0.  It  follows  immediately  now  from  a  theorem 
of  Bernstein  (1929)  and  Widder  (1931)  (cf  Tamarkin  1931)  that 

G(z)  =  Je-*'  dP.it)  (16) 

0 

where  F,(/)  is  a  bounded  nondecieasing  function.  Without  loss  of  generality,  we  may  assume 
^i(0”)  =  0,  (if  not,  sim{dy  add  a  constant,  -F,((r)).  Since  G  (0;  A/)  =  \,N  >  0,  it  follows  immedi¬ 
ately  that  Ff{oo)  =  1.  Finally,  since  is  continuous,  ihe  value  of  F,(r)  at  any  discontinuities,  of 
which  there  can  be  at  most  a  countable  number  (Kolmogorov  and  Fomin  1970),  may  be  chosen  so  as 
to  make  F,(r)  right-continuous  without  affecting  the  value  of  G(z).  Hence,  F,(t)  is  a  CDF  of  a 
positive  random  variable,  and  the  theorem  is  proved. 

2J2  Canonical  Point  Mass  Functions 

The  limit  theorem  given  above  is  predicated  on  the  convergence  of  the  functions  G(z;  N)  to 
G  (z ).  Therefore,  it  is  of  interest  to  study  under  what  conditions  convergence  occurs.  To  investigate 
the  structure  of  the  collections  of  PMF’s,  (pin;  N),N  >0).  that  yield  convergence  in  the  above 
development,  the  following  corollary  to  a  theorem  of  Boas  (1939)  appears  to  be  useful. 

COROLLARY:  Let  be  a  convei^ent  series.  Then  there  exists  a  function 

«  =  0 

p;  (0,  «*>]->/?  of  bounded  variation  such  that 


n  =  0,  1,  2, .  .  . 


(17) 


PROOF:  By  Boas’  theorem,  it  is  permissible  to  write 


a„  =  j^dF(t)  /I  =0,1,2 -  (18) 

0  '*• 


N  jft 

where  F(t)  is  a  function  of  bounded  variation.  Define  //v(t)  =  X  — r»  ^  ®  /n(^) 

n  =  0  ^  ‘ 

increases  with  N  to  e‘  and  by  the  monotone  convergence  theorem 


=Jc'  dFit)  =  C,  (19) 

«  =  0  0 


which  is  finite  by  hypothesis. 
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Define 

t 

m  =  ie^dFia) 

0 

(  t 

=  J  dF^ia)  -  j  e“  dF 2(.a) 

0  0 


=  3i(/)-p2a),  (20) 

where  Flit)  and  F2(/)  are  monotonic  nondecreasing  functions  of  r.  It  follows  immediately  that  3i(0 
and  32(^)  ^  monotonic  nondecreasing  functions  of  t,  and  since  by  Eq.  (19) 

lim  p(f)  =  C  <oo,  (21) 

l-X>o 

it  follows  that  Pi(/)  and  ^2(0  ^rc  bounded:  hence  3(0  is  a  function  of  bounded  variation. 

Now  examine 

t 

3(f)  =  Je“dF(a).  (22) 

0 

where  without  loss  of  generality  we  may  assume  that  F{t)  is  a  monotonic  nondecreasing  function  of 
t.  Also  let  3(f)  =  C,  /  <  0.  Let  ftp  and  ii/;  be  the  Lebesgue-Stieltjes  measures  associated  with  3  and 
F ,  respectively.  Finally,  let  F®  be  the  field  consisting  of  finite  disjoint  unions  of  left-open,  right- 
closed  intervals,  (a,  b],  a  <  b,  a,  b  e  R .  By  suJghtforward  computation,  one  may  easily  show 
that 

lip(>i)  =  rfiif(a)  (23) 

A 

for  all  sets  A  e  F®.  Since  each  side  of  Eq.  (23)  is  a  measure  on  F®,  and  each  of  these  measures  has  a 
unique  extension  to  the  minimal  a-field  over  F®,  which  in  this  case  consists  of  the  Borel  sets,  it  fol¬ 
lows  that  Eq.  (23)  holds  for  all  Borel  sets.  By  the  Radon-Nikodym  theorem, 

li/7  (A )  =  I  e  d  iip(a),  (24) 

A 
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and  hence  from  Eq.  (18) 


=  f  n  =  0,  1,  2,  .  .  .  .  (25) 

6 

This  com[4etes  the  proof  of  the  corollary. 

In  particular,  the  corollary  states  that  any  PMF,  p(n;  N),  n  =  0,  1,  2,  .  .  .  ,  has  the  representa¬ 
tion  given  by  Eq.  (17),  i.e.,  any  random  variable  with  such  a  PMF  is  a  "mixture"  of  Poisson  random 
variables.  Use  of  this  representation  in  the  definition  of  G(z)  and  G(z;  N)  now  yields  in  a  straight¬ 
forward  manner 


G(z)=  lim  J  e  "  dp^(M),  (26) 

0 

where  is  a  function  of  bounded  variation  associated  with  the  representation  of  p(n;  ^).  Little 
may  be  said  about  the  existence  of  G(z)  as  defined  by  Eq.  (26)  without  better  knowledge  of  P’j7(f)- 
However,  it  appears  possible  to  define  a  set  of  parameterized  PMF’s  that  yield  convergence  of  G(z) 
and  that  are  canonical  in  the  sense  that  if  all  collections  of  PMF’s  that  yield  convergence  to  the  same 
function  are  identified  in  an  equivalence  class,  then  this  canonical  set  of  parameterized  PMF’s  contains 
one  representative  of  each  possible  equivalence  class.  In  that  sense,  the  canonical  set  represents  all 
collections  of  PMF’s  that  yield  convergence. 

In  the  theorem,  let  k  =  1  and  let  be  the  variance  of  the  X,  ’s.  The  result  above  becomes 
Cy(u)  =  G (o^u^  /  2).  From  the  relation  =  --Cy^^\0)  is  easily  found  the  relation  G^*\0)  -  -1, 
which  implies 


J  t  dF,(t)  =  1.  (27) 

0 

Thus,  any  CDF  satisfying  Eq.  (27)  yields,  through  Eq.  (16),  an  admissible  G(z)  for  this  problem.  On 
the  other  hand,  if  F,(/)  does  not  satisfy  Eq.  (27),  the  resulting  G(z)  obtained  by  Eq.  (16)  is  not 
admissible  for  this  problem.  For  each  CDF  F,  such  that  F,(-<»)  =  F,(0~)  =  0,  Ft(°»)  =  1,  and  Ff 
satisfies  Eq.  (27),  define  a  collection  of  PMF’s  parameterized  by  N  >  0  by 

oo 

Pc(n;lV)  =  /  -^e-'  dF,(tlN)  n  =  0,  1,  2 .  (28) 

6  n! 

The  canonical  set  PM{N)  (for  Poisson  mixture  with  mean  N)  is  defined  to  be  the  set  of  all  such 
parameterized  PMF’s  as  defined  by  Eq.  (28).  It  follows  from  Eq.  (26)  that  if  ^PM(N),  then 
G(z)  associated  with  p^  is  given  by  Eq.  (16). 
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22.1  Examples 

Let  N  be  a  Poisson  random  variable  with  PMF 

Pin,  N)  =  N>0.  (29) 

n! 

It  is  immediately  evident  that  this  parameterized  PMF  is  a  member  of  the  canonical  set  with 
Ff(t)  =  U it  -  1),  where  Uit)  is  the  Heaviside  step  function.  Hence,  C7(z)  =  e~^,  and  Sfj  con¬ 
verges  to  a  normal  random  variable.  To  exhibit  other  members  of  the  equivalence  class  represented  by 
the  Poisson  distribution,  examine  the  degenerate  case,  with  PMF  given  by 

pin;  N)  =  1,  n  =N  >0,  (30) 


and  the  binomial  distribution,  with  PMF  given  by 

p(/«;^)=  «=  0,1,2 . N,  (31) 

where  0<p  =  l-  ^<l,  andA(  =  Np.  That  the  degenerate  case  is  a  member  of  this  equivalence 
class  is  merely  a  statement  of  the  classical  form  of  the  CLT.  It  is  straightforward  to  show  that 
G  (z)  =  in  this  case.  To  see  that  the  binomial  distribution  is  also  in  this  equivalence  class,  exam¬ 
ine  G  (z ;  N ).  After  a  simple  computation,  one  finds 


Giz;N)  =  il-^f. 

N 

For  a  fixed  p  >  0,N  — » <»  if  and  only  if  ^  — >  <».  It  follows  then  that 


(32) 


limG(z;iV)=  Um  (1  -  -^)^ 
>oo  N 


=  e  \ 


(33) 


which  shows  that  the  binomial  distribution  is  also  in  this  class. 

As  a  second  example,  let  N  be  a  negative  binomial  random  variable  with  PMF 


/  r(n-»-v)  - 


(34) 


where  r  =  N/v.  This  PMF  is  also  a  member  of  the  canonical  set  with 
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=  (35) 

0  * 

as  may  be  verified  directly.  Hence  G(z)  converges  for  the  negative  binomial  distribution  and  is  given 
by 


G(z)  = 


V 

v+z 


(36) 


The  attributes  of  the  canonical  PMF’s  that  are  of  interest  in  this  problem  are;  a)  they  forni  a  set 
of  parameterized  PMF’s  for  which  convergence  always  occurs,  b)  they  represent  all  collections  of 
PMF’s  for  which  convergence  occurs  and  c)  they  are  related  to  the  limiting  CDF’s  F,(f )  in  a  straight¬ 
forward  way. 

3.  CLUITER  MODELLING:  RAYLEIGH  MIXTURE  DISTRIBUTIONS 

If  the  random  vectors  in  the  theorem  are  defined  in  R^,  the  limiting  random  variable  denoted  as 
F  =  (y  1  Y 2)  and  the  matrix  Q  =  6^1  where  >  0  is  a  constant  and  /  is  the  2  x  2  identity  matrix, 
then  the  problem  being  discussed  here  takes  the  form  of  a  spherically  symmetric  random  walk  in  the 
plane  with  a  random  number  of  steps.  Typically  in  ai^lications  involving  narrowband  signals  (or 
noise),  the  two  components  Y  ^  and  Y 2  iic  interpreted  as  the  inphase  and  quadrature  components  of  a 
complex  random  variable.  The  relationship  with  the  phasor  approach  used  in  formulating  Eqs.  (1)  and 
(2)  is  established  by  transforming  Yi  and  Y2to  polar  coordinates,  say  R  and  <{>.  The  assumption  that 
is  tantamount  in  this  problem  to  assuming  that  R  and  (]>  are  independent  and  that  (|>  is  uni¬ 
formly  distributed  on  [0,  27c].  If  the  resultant  (or  amplitude),  /?,  of  the  limiting  case  of  this  random 
walk  isN  — >  00  is  defined  by  /?^  =  Y^^  +  Y‘^,  then  inversion  of  the  characteristic  function  yields 


/«<*)  =  ^  dF,(t)  (37) 

where  Flit)  is  defined  in  the  theorem.  The  choice  of  F,(r)  obtained  from  the  negative  binomial  dis¬ 
tribution  yields  the  result  given  by  Jakeman  and  Pusey  (1978).  Evidently,  the  only  permissible  distri¬ 
butions  for  this  resultant  are  mixtures  of  Rayleigh  distributions.  Use  of  this  type  of  model  for  the 
amplitude  statistics  of  clutter  has  been  proposed  by  previous  investigators,  whose  work  is  discussed 
below.  The  derivation  of  this  model  that  is  presented  here,  however,  is  new. 

After  studying  the  problem  of  sea  clutter  modelling.  Trunk  (1972;  1976)  concluded  that  some 
types  of  non-Rayleigh  sea  clutter  may  be  modelled  as  a  locally  homogeneous  Rayleigh  process  whose 
Rayleigh  parameter  (which  represents  clutter  power  in  this  case)  is  modulated  by  the  radar’s  large  scale 
spatial  sampling  of  the  environment.  The  amplitude  PDF  of  such  a  model  is  given  by  Eq.  (37).  In  the 
remainder  of  this  report,  the  model  described  by  Eq.  (37)  is  referred  to  as  a  Rayleigh  mixture  model. 
Prior  to  this  work,  Valenzuela  and  Laing  (1971)  had  argued  that  the  scattering  of  a  radar  pulse  from 
the  surface  of  the  sea  modelled  by  the  so-called  composite  surface  model  yields  amplitude  statistics 
whose  PDF,  although  complicated  in  form,  may  be  described  as  a  Rayleigh  mixture  model.  Jakeman 
and  Pusey  (1976)  used  the  phenomenological  model  described  by  Eq.  (1)  to  study  radar  sea  clutter 
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modelling.  Specifically,  they  assumed  the  ‘s  to  be  i.i.d.  according  to  a  K  distribution,  the  <]>,■  ’s  to  be 
i.i.d.  according  to  the  uniform  distribution  on  [0.  2k],  and  the  <2,-  ’s  to  be  mutually  independent  of  the 
<!>,■  ’s.  FitMn  these  assumptions,  they  showed  that  die  amplitude  of  the  scattered  field  E  is  statistically 
described  by  a  K  distribution,  which  is  also  an  example  of  a  Rayleigh  mixture  model.  Note  that  in  the 
1976  work,  the  number  of  scatterers  N  was  a  deterministic,  nonfluctuating  quantity;  hence,  the  need 
for  a  ^cific  model  of  the  statistics  of  the  individual  scatterers  a,-  arose.  In  the  1978  work  cited 
above,  they  presented  an  alternative  derivation  of  K  distributed  amplitude  statistics  by  introducing 
number  fluctuations  into  the  model  and  relaxing  the  requirement  that  the  individual  's  have  K  disui- 
buted  statistics.  Since  the  introduction  of  the  K  di^bution,  it  has  been  used  with  qualitative  success 
in  empirical  fits  to  clutter  data  (Ward  1981,  1982). 

In  addition  to  the  models  described  above,  the  WeibuU  distribution  has  also  yielded  qualitatively 
accurate  empirical  fits  to  clutter  data  (Skolnik  1980;  Olin  1982;  Trizna  1988,  1989)  and  has  emerged  as 
a  candidate  model  for  sea  clutter  statistics.  Its  effectiveness  as  a  model  appears  to  have  been  strictly 
empirical;  it  does  not  appear  to  have  the  physical  appeal  of  the  K  distribution  or  other  Rayleigh  mix¬ 
ture  models.  Conte  and  Lx)ngo  (1987)  addressed  this  issue  through  an  analysis  that  strongly  suggests 
that  the  Weibull  distribution  is  also  a  Rayleigh  mixture  distribution  and,  hence,  would  have  the  same 
physical  interpretation  as  the  K  or  other  Rayleigh  mixture  models.  Unfortunately,  they  were  not  able 
to  reach  a  definitive  conclusion.  A  demonstration  based  on  an  early  result  of  Bochner  (1937a)  is 
presented  below  and  shows  rigorously  that  the  Weibull  distribution  is  in  fact  a  Rayleigh  mixture  distri¬ 
bution  for  a  limited  range  of  the  shape  parameter  (i.e.,  b  <  2;  see  Eq.  (47).) 

Thus,  as  hrst  suggested  by  Conte  and  Longo,  almost  all  of  the  major  statistical  models  that  have 
been  proposed  for  the  amplitude  statistics  of  sea  clutter  (the  log-normal  being  an  exception)  may  be 
represented  as  a  Rayleigh  mixture  model.  This  observation  coupled  with  the  result  that 
physical/phenomenological  models  that  have  been  proposed  to  model  the  scattering  process  also  lead  to 
the  Rayleigh  mixture  model  suggests  that  the  model  described  by  Eq.  (37)  be  considered  a  unified 
model  for  sea  clutter  amplitude  statistics.  A  study  of  the  properties  of  such  distributions,  especially  to 
determine  if  any  of  the  previously  proposed  models  for  amplitude  statistics  of  scattering  problems  are 
examples  of  Rayleigh  mixture  distributions,  is  therefore  warranted  and  is  given  below. 

3.1  Conditions  for  Rayleigh  Mixture  Distributions 

A  CDF  defined  on  the  positive  real  line  will  be  said  to  be  the  distribution  of  a  Rayleigh  mixture 
if  this  CDF,  {x ),  may  be  represented  in  the  form 


F,(x)  =  /(!-€  2^)dF,(x).  (38) 

0 

where  x  (i.e.,  the  Rayleigh  parameter)  is  a  random  variable  defined  on  the  positive  real  line,  and  F^(x) 
is  the  CDF  of  X.  The  problem  of  interest  here  is  to  determine  necessary  and  sufficient  conditions  on 
FjfOc)  to  ensure  that  it  has  the  representation  given  in  Eq.  (38)  for  some  CDF,  F^(x).  These  condi¬ 
tions  follow  as  a  corollary  of  the  Bernstein- Widder  theorem  for  c.m.  functions  that  was  referred  to  ear¬ 
lier  in  the  proof  of  the  theorem  presented  herein. 

COROLLARY:  A  CDF  F^ix)  defined  on  the  positive  real  line  is  a  Rayleigh  mixture  distribution 
if  and  only  if  the  associated  function  K(y)  =  1  -  F,  (Vy  )  is  c.m.  in  the  interval  0  <  y  <  «». 
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PROOF:  Assume  i<r(y)  is  c.m.  in  the  interval  0  <  «»  .  Then  by  the  Bernstein- Widder 

theorem. 


K(y)  =  je-y^dP(a),  (39) 

0 

where  ^(a)  is  ncxidecreasing.  Without  loss  of  generality.  ^(0)  may  be  chosen  to  be  equal  to  0.  This 
choice  has  no  effea  on  the  value  of  K(y).  Since  ^^(0)  =  1,  it  follows  that  P(<»)  -  P(0)  =  pC**)  =  1. 
P(a)  is  seen  to  be  bounded  and,  therefore,  has  at  most  countably  many  points  of  discontinuity.  Since 
e~y°'  is  continuous,  the  value  of  P(a)  at  each  point  of  discontinuity  may  be  assigned  to  be  an  aibitrary 
finite  value  without  affecting  the  value  of  K(y).  Therefore,  the  values  of  P(a)  at  its  points  of  discon¬ 
tinuity  can  be  chosen  so  as  to  make  P(a)  right-continuous. 

Clearly,  y  =x^  yields 


F,(x)  =  I  (1  -  c-"**)  dP(a). 


(40) 


The  substitutions  a  =  l/2t  and  F^(x)  =  1  —  P(l/2x)  yield 


0 


(41) 


It  remains  only  to  show  diat  F^x)  is  a  CDF.  To  this  end. 


1.  F^(x)  is  nondecreasing,  since  P(a)  is  nondecreasing. 

2.  lim  F^x)  =  Urn  F^x)  =  0  since  lim  P(a)  =  1. 

T-*  -o»  t-*o  o-»«~ 

Similariy,  lim  F,(x)  =  1. 

T-»  «■ 

3.  F.{(x)  is  right-continuous,  since  P(a)  is  right-continuous. 


Therefore,  F^ix)  is  a  CDF,  and  the  proof  in  one  direction  is  complete. 

Conveisely,  assume  Fj^(x)  is  a  Rayleigh  mixture  distribution.  Then,  from  the  substitutions 

a  =  l/2x  (42) 


and 


p(a)  =  1  -  F,(l/2a). 


(43) 
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it  foUows  that 


K{y)  =  \e-y^d^ia)  (44) 

0 

where  P(a)  is  nondecreasing.  Therefore,  by  the  Bemstein-Widder  theorem,  /iT  (y )  is  c.m.  in  the  inter¬ 
val  0  <  y  <  o®.  Because  #r(0)  =  K{Q^)  =  1,  K(y)  is  also  c.m.  in  the  interval  0  ^  y  <  <».  The 
proof  is  complete. 

This  corollary  fully  characterizes  Rayleigh  mixture  distributions  and  establi^es  a  one-to-one  rela¬ 
tionship  between  Rayleigh  mixture  CDF’s  and  functions,  K{y),  that  arc  c.m.  in  the  interval 
0  ^  y  <  oo  and  normalized  to  ^^(0)  =  1.  Given  a  CDF  one  may  establish  that  it  represents  a 
Rayleigh  mixture  distribution  by  establishing  that  the  associated  function  AT  is  c.m..  Although  this  task 
requires  examination  of  an  infinite  number  of  derivatives,  some  techniques  are  presented  below  that 
make  this  problem  tractable  in  some  instances.  Further  conditions  that  ensure  dF^(x)  =  /x(t)  dx  fol¬ 
low  from  the  work  of  Widder  (1931)  and  Hille  and  Tamarkin  (1933,  1934).  The  interested  reader  is 
referred  to  the  references  for  details. 

32  Examples  of  Rayleigh  Mixture  Distributions 

To  assess  if  a  given  distribution  is  a  Rayleigh  mixture,  two  previously  obtained  results  relating  to 
the  class  of  c.m.  functions  are  useful.  First,  let  C  be  the  class  of  functions  that  are  c.m.  on 
0  <  y  <  oo,  Bochner  (1937a)  shows  that 


1-  If  / 1»  /2,  E  ^  •  ihen  a/ 1  +  4/2  e  C  fora,  b  >  0. 
2.  If  / 1,  /2,  £  C 1  then / 1/2  £ 


The  second  result  is  a  lemma  also  due  to  Bochner  (1937a,  pp.  498-499); 

LEMMA:  If  \|t(p)  is  c.m.  on  0  ^  p  <  00  and  x(p)  is  such  that  x(0)  =  0,  and  X^*^(P)  is  c.m.  on 
0  <  p  <  00,  then  the  composition  \|t(X(p))  is  c.m.  on  0  S  p  <  o®. 

As  stated  earlier,  a  well-known  example  of  a  Rayleigh  mixture  distribution  is  the  K  distribution 
with  PDF  defined  by 


fx(x)=  ^2^1  r(v)  ^ ^  X ,  V,  Ti  >  0, 


from  which  the  following  representation  may  be  obtained  (Watson  1944,  p.l83) 


r(v) 


-  — t 


(45) 


(46) 
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A  second  examine  is  given  by  the  Weibull  distribution  with  PDF 

f^ix)  =  abx^~^e~*“*  x,a,b>0.  (47) 

If  the  parameter  b  is  restricted  to  the  range  0  <  b  ^2,  Bochner’s  lemma  may  be  applied  to  show  that 
this  PDF  is  a  mixture  of  Rayleigh  PDF’s  (Bochner  1937b).  To  see  this  result,  examine 

K(y)=l-GC^) 


=  e-^  .  (48) 

In  Bochner’s  lemma,  let  V(p)  =  e~^^,  which  is  easily  seen  to  be  c.m.  on  0  ^  p  <  oo  for  any  a  >  0, 
and  let  x(y)  =  Clearly,  x(0)  =  0  and  =  ib/2)y^^~^,  which  for  <2  may  be  written 


^tn-i  ^ 
2^ 


bl2 

ni  -bn.) 


J  dt. 
0 


(49) 


By  the  Bernstein- Widder  theorem,  Eq.  (49)  shows  that  (bn)y^^~^  is  c.m.  on  0<}'  <«»  for 
0<b  <2.  By  Bochner’s  lemma,  then,  ^r(y )  is  c.m.  on  0  y  <  «».  it  follows  that  the  Weibull  dis¬ 
tribution  is  a  Rayleigh  mixture  distribution  for  the  range  0  <  b  !^2.  (The  case  ^  =  2  is  obvious.) 

To  tqjply  the  results  presented  in  this  repott  to  a  lesser  known  case,  consider  as  a  third  example 
the  gamma  distribution  with  PDF 


fx(x)  = 


nb) 


X,  a,  b  >  0. 


(50) 


Oeariy,  the  following  relations  hold: 


K(y)=l-F^(^)^0, 


(51) 


_/j:(l)(y)  =  L01  ^  0.  (52) 

Also,  since 

(_l)«A:(«)(y)  = 

dy"~^ 

2Vy 
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-1 


dy 


n—l 


(53) 
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to  show  ^  0,  it  suffices  to  show 


(-1) 


,n-l 


rt-1 


dy 


n-1 


2<y 


^  0  n  ^  1, 


(54) 


i.e.,  that  the  function  )  /  2^7  is  c.m.  on  0  <  y  <  <».  To  show  the  complete  monotonicity  of 
this  function,  examine 


2^7  br(b)  2  ^ 


(55) 


It  suffices  to  show  that  the  functions  and  (bl2)y^'^~^  are  each  c.m.  and  to  use  the  closure  pro¬ 

perties  of  c.m.  functions.  But  these  fim.tions  have  already  been  shown  above  to  be  c.m.  for 
0  <  b  <2.  If  b  =2,  then  (bl2)y^'^~^  =  1,  which  is  c.m.  because  all  of  its  derivatives  are  equal  to  0 
(indeed,  any  positive  constant  is  c.m.  for  the  same  reason).  Therefore,  by  the  closure  property  of  c.m. 
functions,  the  gamma  PDF  is  also  a  mixture  of  Rayleigh  PDF’s  for  the  case  0  <  b  <2. 

As  a  final  example,  consider  the  Nakagami-m  distribution  with  PDF  given  by 


2m-l  - 

/x(^)  =  2(— r  -  ■■■  e  ,  x,m,r\>0.  (56) 

T\  r(m) 

A  computation  similar  to  the  one  above  for  the  gamma  distribution  shows  that  this  distribution  is  a 
Rayleigh  mixture  distribution  for  m  <1. 

33  Exponential  Mixture  Distributions 

As  it  is  sometimes  of  interest  to  study  the  distribution  of  intensity  rather  than  the  distribution  of 
amplitude,  let  /  =  /?^  in  Eq.  (37).  In  radar  applications,  I  represents  the  radar  cross-section.  The 
CDF  of  intensity  becomes 


_  JL 

/^/(>’)  =  /(l-e  2")dF,(T),  (57) 

0 

where  t  =  ah  and  F-^{x)  =  F,  (z/ah-  Distributions  of  this  form  will  be  called  exponential  mixture 
distributions. 

In  problems  of  propagation  through  random  media,  previous  investigators  have  sometimes  studied 
the  moments  of  intensity  rather  than  its  distribution.  From  Eq.  (57),  the  moments  of  intensity  predicted 
by  this  model  are  given  by 


=2''n!Elz^]. 

</>" 


(58) 


This  model  is  consistent  with  models  presented  by  previous  investigators  (Dashen  1984).  From  Eq. 
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(58),  (xie  inunediately  obtains  for  the  scintillation  index  S , 

</2> 

S  =  -1^1.  (59) 

<I>^ 

with  equality  if  and  only  if  F^(x)  is  a  unit  step  ftmction.  For  application  to  problems  of  propagation 
through  continuous  random  media,  this  model  is  thus  only  applicable  in  the  saturation  regime.  But  this 
region  is  precisely  the  region  in  which  Dashen  (1979)  indicates  that  the  phasor  interpretation  described 
by  Eq.  (2)  is  ^plicable. 

4.  METHODOLOGIES  FOR  DETERMINING  F,(t) 

4.1  Theoretical  Inversion  of  a  Rayieigh  Mixture  Distribution 

Thus  far  the  results  concerning  Rayleigh  mixtures  have  ascertained  only  the  existence  of  the  CDF 
Once  the  Rayleigh  mixture  property  is  established  for  a  given  f^ix),  the  natural  question  to 
ask  is  what  is  the  form  of  the  associated  F^(x)7  The  general  problem  appears  to  be  quite  difficult.  As 
a  result,  the  discussion  in  this  section  is  limited  to  suggestions  for  attacking  the  problem. 

A  general  result  for  the  inversion  of  Eq.  (39),  i.e.,  obtaining  P(a)  for  a  given  AT  (y ),  is  presented 
by  Widder  (1934).  Note  that  F,(t)  may  be  obtained  directly  from  P(a)  by  a  simple  change  of  vari¬ 
able.  For  the  problem  of  interest  here,  the  result  becomes 

p(a)  =  lim (-!)«  f  du.  (60) 

n-*oo  i  n! 

n 

a 

Although  general,  the  result  requires  knowledge  of  the  derivatives  of  all  orders  of  AT  (y )  in  a  neighbor¬ 
hood  as  y  — >  oo  .  Since  the  actual  computation  of  these  derivatives  can  be  difficult  for  some  choices 
of  Ar(y ) ,  this  method  can  be  intractable.  The  interested  reader  is  referred  to  the  reference  for  details. 

A  second,  quite  general  result  for  the  inversion  of  Eq.  (39)  is  also  given  by  Widder  (1935).  This 
method  requires  knowledge  of  the  derivatives  of  all  orders  of  AT  (y )  at  the  point  y  =  1  and  yields  the 
following  series  solution  for  P(a)  in  terms  of  Laguerre  polynomials: 


This  approach  may  be  illustrated  by  an  application  to  determine  the  inversion  for  the  WeibuU  distribu¬ 
tion.  In  this  case,  K(y)  =  where  c  =  b/2.  To  apply  Widder’s  result,  knowledge  of  AT^^^l)  is 
necessary.  These  derivatives  do  not  appear  to  be  obtainable  in  closed  form.  However, 


dfy""* 


(62) 
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Application  of  the  Leibniz  nile  for  the  differentiation  of  a  product  then  yields 

^(»»)(y)  =  -acyc-n  2)  1^”  "  n  ^  1, 


where  c*  =  (c  -  1)  (c  -  2)  •  •  •  (c  -  n  +  ik  +  1).  This  computation  is  a  recursive  computation 
for  the  derivatives  of  K(y).  It  therefore  may  be  used  in  conjunction  with  Widder’s  result  to  give  a 
recursive  scheme  for  computing  P(a)  and,  thus,  The  drawback,  of  course,  is  that  this  reproach 

does  not  result  in  a  closed  form  solution  for  F,(x). 

For  many  cases  of  interest,  X  will  be  a  continuous  random  variable  and  dF^(x)  =  /^Ct)  dx  will 
hold.  In  cases  for  which  this  relationship  does  hold,  Eq.  (39)  is  an  ordinary  Laplace  transform,  and 
/.j(x)  may  be  obtained  (at  least  conceptually)  from  the  inverse  Laplace  transform  of  Eq.  (39).  This 
inversion  result  is  given  by 

1  1  °  ^  )  — 

However,  in  some  cases  inversion  via  the  standard  contour  integral  approach  may  be  difficult  For 
example,  for  the  Weibull  distribution,  Eq.  (65)  becomes 

/tW  =  ds  a  >  0,  (66) 

where  c  =  b/2.  The  multivaluedness  of  the  integrand  renders  this  approach  quite  difficult  (see,  how¬ 
ever,  PoUard  (1946)). 

Under  the  assumption  that  /.^(x)  is  piecewise  continuous,  an  alternate  approach  for  computing 
fxi^)  based  on  the  Mellin  transform  (Courant  and  Hilbert  1953)  is  as  follows. 

Consider  the  function 


0 


If  the  integral  converges  for  Re(p)>  C,  then 


o® 


+  1)  =  J  I  -c  2v,(x)  dxdx  Re(p)  >  o. 


0  0 
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The  integrals  on  the  right-hand  side  may  be  interchanged  to  obtain 


F^j^Clp  1)  =  2PT(p  +  1)/  x^/,(x)  dx 


=  2Pr(p  -I- l)F^(p  -t- 1) 


(69) 


where  F^(p)  is  the  Mellin  transform  of  /^('C)  and  is  convergent  for  Re(p)  >  a.  Since  /^(x)  is 
piecewise  smooth  by  assumption,  it  follows  that 


<j  +  j‘ 


/,(t)  =  f  +  Dt-'-'  dp. 


Substitution  for  F^(p  +  1)  yields 


2kj„Jj^  2Pr(p  +  l)  ^ 


(70) 


(71) 


This  result  indicates  that  under  the  appropriate  conditions  (i.e.,  the  piecewise  continuity  of  /,(x)  and 
the  convergence  of  (2p+l)),  /,(x)  as  dehned  by  Eq.  (71)  is  an  inversion  integral  for  this  prob¬ 
lem,  This  inversion  integral  may  also  be  intractable,  but  it  does  provide  an  alternative  approach  to  the 
inversion  integral  of  Eq.  (6S). 

For  an  2q)plication  of  this  approach,  let/jj(x)  be  the  PDF  of  the  K  distribution.  For  this  case,  X 
is  known  to  be  a  continuous  random  variable,  and 


F,J^f(2p +!)  =  ( 


2  .2nr(v  +  p)r(p  +  1) 

V2v/r[  r(v) 


(72) 


This  integral  converges  for  Re(p)>  (l/2)(l/?«(v  -  1)1  -  V  -  1).  Inserting  this  expression  into  Eq. 
(71)  and  letting  z  =  v  +  p  yields 


/.(t)  = 


(V/Tl)^  ,v-l 

r(v) 


a  +  ;« 

—  f  r(z)(vx/Ti)-*  dz 


a>(l/2)(l/?e(v- l)l-»-v- 1).  (73) 


Note  that  the  contour  lies  in  the  right  half-plane  for  any  choice  of  v  >  0.  The  singularities  of  the 
integrand  are  simple  poles  due  to  the  presence  of  r(z )  and  occur  at  z  =  -n ,  n  =  0,  1,  2,  .  .  .  .  The 
residues  of  the  gamma  function  at  these  points  are  equal  to  (-I)"//! !.  Closing  the  contour  to  the  left 
by  a  semicircle  of  radius  Rff  that  encloses  N  poles  and  taking  the  limit  as  — >  <»  yields 
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Vvv 


/t(x) 


r(v) 


pV-1 


z 

i  =0 


_ !1__ 

/! 


(V/Tp^ 

r(v) 


(74) 


which  is  known  to  be  correct. 

An  application  of  this  same  method  to  the  case  in  which  (x )  is  the  PDF  of  the  Weibull  distri¬ 
bution  yields 


F,j^(2p  +  l)  = 


(75) 


The  computaticMi  proceeds  similarly  to  the  one  just  given  and  yields 


f^(x)  =  ab(2x)^ 


£  (-a(>/^)^y 

‘=o,!r(l-|-(i  +1)) 


0<b  <2. 


(76) 


This  result  agrees  with  a  result  obtained  independently  by  Pollard  (1946)  and  by  Grosjean  et  al.  (1989), 
and.  in  contrast  to  the  result  obtained  earlier  herein,  is  a  closed-form  solution.  (The  issue  as  to 
whether  this  result  is  more  desirable  from  a  computational  point  of  view  than  the  previous  result  is  not 
addressed  here.) 

Formal  application  of  this  approach  to  the  case  in  which  (jc )  is  the  gamma  PDF  yields 


f  (2p  +  i)  =  iee-±^. 

a^nb) 


(77) 


Again,  the  computation  proceeds  similarly  to  the  one  given  above  and  yields 


(-aV2T)* 


nb) 


1 


0<b  ^2. 


(78) 


‘  =o,!r(l  -  -^(j  +b)) 


For  b  =  1  and  b  =  2,  this  series  sums  to 


/x(x)  = 


a 

VnVIx 


e 


21 


(79) 


SANGSTON  AND  GERLACH 


and 

respectively  (Prudnikov  et  al.  1986,  p.  707).  The  results  in  Eqs.  (79)  and  (80)  may  be  ascertained  to 
be  correct  inversions  by  direct  computation  via  Eq.  (65),  or  by  inserting  them  in  Eq.  (37)  with  T  =  ah 
and  performing  the  integration. 

Finally,  formal  application  of  this  approach  to  the  Nakagami-m  distribution  proceeds  analogously 
to  the  above  derivations  and  yields 


/,(T)  =  2(^r 


“  i-2xm/r\y 
Tim)  +0)’ 


m  <  1. 


(81) 


42  Fitting  a  Rayleigh  Mixture  to  Data 

If  the  Rayleigh  mixture  model  is  to  used  to  model  data,  the  problem  of  empirically  fitting  data  to 
a  distribution  becomes  a  problem  of  finding  the  F^{x)  that  yields  the  closest  fit  to  the  data.  Moreover, 
since  in  this  case  an  entire  function  rather  than  a  limited  number  of  parameters  (usuaUy  two)  is  chosen 
during  the  fitting  process,  this  model  should  give  closer  empirical  fits  to  the  data.  For  this  approach  to 
be  effective,  however,  a  methodology  for  obtaining  F,(x)  from  a  given  set  of  data  is  required.  This 
section  of  the  report  presents  such  a  methodology.  If  an  estimate  for  F^{x)  can  be  obtained  from  a 
given  set  of  data,  then  the  PDF  of  the  clutter  amplitude  can  be  estimated  from  Eq.  (37)  with  t  =  ah. 
For  f  n(x)  described  by  Eq.  (37)  it  can  be  shown  that 


I  x'^fnix)  dx  =  2PT{p  +  \)jxP  dF^(x),  (82) 

0  0 

where  p  =  a  +  j(0  and  Fip )  is  the  gamma  function,  provided  the  left-hand  side  exists.  For  the  prob¬ 
lem  of  interest  here,  let p  -  n.  Then  Eq.  (82)  yields 


V^x.2n  =2"^! 


(83) 


where 


iig  2n  =  2nth  moment  of  x. 


=  n  th  naoment  of  x. 


As  the  moments  of  the  amplitude  x  can  be  determined  directly  from  experimental  data,  the  moments  of 
t  can  be  calculated  from  Eq.  (83). 
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To  illustrate  this  idea,  consider  in  place  of  real  data  a  model  distribution  for  which  the  moments 
can  be  written  in  closed  form.  For  (x )  modelled  by  the  K  distribution,  the  moments  are  given  by 


=  (^y 


r(v) 


which  yields  for  the  moments  of  X 


r(v  +  n) 
r(v) 


(84) 


(85) 


These  moments  are  the  moments  of  the  gamma  distribution,  which  is  known  to  be  the  correct  Fj^(x) 
for  g(x)  given  by  the  K  distribution. 

The  problem,  therefore,  is  to  find  a  distributitm  F^ix)  that  is  consilient  with  a  set  of  moments 

obtained  as  described  above.  This  problem  is  of  classical  origin,  and  many  approaches  are  avail¬ 
able.  One  approach,  based  on  an  orthogonal  polynomial  expansion,  will  be  examined  here.  In  this 
approach,  it  is  assumed  that  the  distributicn  of  x  admits  a  probability  density  function  /  ^(x).  Although 
this  assumption  seems  rather  arbitrary,  the  empirical  results  obtained  in  the  past,  namely  the  success  of 
the  K  and  WeibuU  distributions,  suggest  that  this  assumption  is  valid  for  modelling  sea  clutter.  With 
this  assumption,  Eq.  (37)  (with  x  =  a^r)  becomes 


^‘'f^x)dx,  (86) 

6  ^ 

and  the  task  is  to  determine  /x(x)  from  the  moments  n  =  0,  ,  N. 

To  determine  this  PDF,  it  is  now  assumed  that/.t(x)  may  be  expanded  in  an  infinite  series  of 
Laguerre  polynomials 


/(t)  =  £  a>-hb>0,  (87) 

« =0  *  ® 

where  L^^ix)  is  the  generalized  Laguerre  polynomial  of  order  n  and  a  and  b  are  constants  to  be 
determined  later.  Since  these  polynomials  are  orthogonal  on  the  interval  [0, «»)  (which  is  the  interval 
of  interest  for  this  problem)  with  respect  to  the  weight  x^e~^,  they  are  an  attractive  choice  as  basis 
functions  for  the  expansion  of  the  uttknown  f^ix).  This  orthogonality  condition  implies 


= 


^  • _ f  /  a(2.)£S^  dx 

r(a  +  H-n)i  6 


(88) 
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upon  substimtion  of  the  foimub  for  the  Lagueire  polynomial,  becomes 

n\  "  (-ir  fn  +  a] 

"  r(a  +  1 +  n)^"o  U-rnl^m+r 


(89) 


The  ooefhcients  in  the  expansion  are  seen  from  Eq.  (89)  to  be  directly  obtainaUe  from  knowledge  of 
the  moments  of  x.  Once  this  expansion  has  been  obtained,  it  may  be  substituted  into  Eq.  (86)  to 
obtain 


/,(x)  =  2i  £  c.  i  (90) 

n=0  m=0  ml  VI  rrv  2 

The  parameters  a  and  b  are  arbitrary  parameters  that  define  the  Laguerre  polynomials  used  in  the 
expansion  and  may  be  advantageously  chosen  as  follows.  Knowledge  of  the  moments  up  through  the 
second  order  implies  that  the  expanaon  has  three  terms  widi  the  coefficients  Cq  through  C2  depending 
on  these  three  moments  (note  that  |io  =  1  is  considered  a  moment).  Because  a  and  b  are  aibitrary,  it 
is  possible  to  choose  these  parameters  so  that  c  ]  =  C2  =  0.  This  choice  is 

a  =  — j  ~  1,  (91) 

m,2  -  m.1 


(92) 


Note  that  in  Eq.  (91)  the  denominator  is  the  variance  of  t,  whereas  in  Eq.  (92)  the  numerator  is  this 
variance.  Denote  this  variance  by  and  let 


(93) 


and 


With  these  choices,  the  expression  for  the  expansion  coefficients  becomes 


r(v  +  m) 


m,w 


(94) 


(95) 
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42.1  Examples 

As  real  clutter  data  was  not  available  for  this  study,  the  application  of  the  methodology  developed 
herein  will  be  illustrated  with  models  for  which  closed  form  expressions  for  the  moments  are  available. 
First  consider  clutter  for  which  the  amplimde  is  K  distributed.  In  this  case,  the  moments  of  T  are 
given  by  Eq.  (85),  and  substitution  in  Eq.  (95)  yields 


(^) 


=  — ^  E  (-ir  • 

nv)^to 


(96) 


The  sum  in  Eq.  (96)  equals  1  for  n  =  0  and  equals  0  for  n  >  0.  Thus,  based  on  knowledge  of 
moments  of  the  K  distribution  up  through  at  least  second  order,  the  approach  described  herein  yields 

//?  U )  =  X ),  (97) 

which  is  the  correct  PDF  for  the  K  distribution. 

Consider  now  a  case  in  which  the  clutter  is  not  exactly  K  distributed  but  instead  is  distributed 
according  to  the  WeibuU  distribution.  Specifically,  for  the  PDF 

g{x)  =  abx^~^e~"\  a,b  >0, 
the  related  moments  of  X  for  use  in  Eq.  (95)  are 

+  1)  . 

"  Tin  +  1)  (2a2/fr)«  ’ 

Figures  1  through  3  compare  the  exact  distributions  with  the  distributions  obtained  from  the  expansion 
method  described  herein.  In  all  cases,  the  series  approximation  converges  to  the  actual  distribution 
function  with  only  a  small  number  of  terms. 

As  a  second  example,  consider  the  case  in  which  the  true  distribution  of  the  clutter  is  a  mixture 
of  Weibuli  distributions  of  the  form 


(98) 


(99) 


G(x)  =  (l  -p)W(x;a,,  Z7,)  +  pW(jc;a2,  ^2).  (100) 


where  W{x;  a,  b)  represents  a  Weibuli  distribution  with  parameters  a  and  For  concreteness,  let 
a  I  =02=  I,  b  I  =  \,  b2  =  2,  and  p  =  0.7.  Physically,  this  distribution  corresponds  to  Rayleigh 
clutter  ib  =  2)  that  is  corrupted  by  spikes  {b  =  1)  30%  of  the  time.  Figure  4  compares  the  exact  mix¬ 
ture  distribution  with  the  distributions  obtained  by  the  expansion  method. 
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0.0  1.0  2.0  3.0  4.0  5.0  6.0 

THRESHOLD 


Fig.  1(a)  —  Approximation  of  WeibuU  (a  =  I,  b  =  1.9)  with  2  moments 
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Fig.  1(b)  —  Approximation  of  Weibull  (a  =  1,  ^>  =  1.9)  with  3  moments 
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Fig.  1(c)  —  Approximation  of  WeibuU  {a  =  I,  b  =  1.9)  with  4  moments 


28 


6.0 


NRL  REPORT  9367 


Fig.  2(a)  —  Approximation  of  WeibuU  (a  =  1,  =  1.5)  with  2  moments 
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Fig.  2(b)  —  Approximation  of  WeibuU  (a  =  1,  ft  =  1.5)  with  3  moments 
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Fig.  2(c)  —  Approximation  of  WeibuU  (a  =  1,  =  1.5)  with  4  moments 
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0.0  1.0  2.0  3.0  4.0  5.0  6.0 

THRESHOLD 

Rg.  2(d)  —  Approximation  of  Weibull  (a  =  1,  ft  =  1.5)  with  5  moments 
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^  I  I  I  I - 1 - 1 

0.0  1.0  2.0  3.0  4.0  5.0  6.0 


THRESHOLD 

Rg.  3<a)  —  Approximation  of  WdbuU  (a  =  1,  6  =  0.8)  with  2  momentf 
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3.0 


4.0 


THRESHOLD 


T“ 

5.0 


Fig.  3(c)  —  Approximation  of  WeibuU  (a  =  \,  b  =  0.8)  with  4  moments 
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THRESHOLD 


Fig.  3(d)  —  Approximation  of  Weibull  (a  =  \,  b  =0.8)  with  5  moments 
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Fig.  4(a)  —  Approximation  of  Weibull  mixture  with  2  moments 
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Fig.  4(h)  —  Approximation  of  WeibuU  mixture  with  5  moments 
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THRESHOLD 


Fig.  4(c)  —  Approximation  of  WeibuU  mixture  with  9  moments 


SANGSTON  AND  GERLACH 


5.  MULTIDIMENSIONAL  NON-GAUSSIAN  CLUTTER  MODELS 

Conte  and  Longo  (1987)  present  the  idea  of  modelling  the  radar  sea  clutter  process  as  a  spheri¬ 
cally  invariant  random  process.  This  idea  is  attractive  as  it  leads  to  a  multidimoisional  model  of  the 
clutter,  namely  a  multidimensional  PDF,  that  includes  correlation  and  need  not  be  Gaussian.  More¬ 
over,  they  also  show  that  most  of  the  major  models  for  non-Rayleigh  amplitude  statisties  may  be 
included  in  the  framework  of  this  model.  As  was  discussed  earlier  in  this  repon,  the  Weibull  distribu¬ 
tion,  for  which  Conte  and  Longo  only  gave  a  partial  result,  is  also  consistent  with  this  framework. 
This  model,  therefore,  has  merit  as  a  multidimensional  model  of  non-Gaussian  clutter  and  deserves 
further  investigation.  In  this  section,  application  of  the  theorem  presented  in  Section  2  is  shown  to 
lead  to  the  same  multidimensional  PDF  that  is  suggested  by  Conte  and  Longo.  Since  the  theorem  ean 
be  interpreted  in  a  physical  sense,  this  result  gives  a  physical  motivation  for  this  approach.  In  addition, 
the  shortcoming  of  the  model  in  terms  of  the  physical  modelling  of  the  problem  emerges  clearly  from 
applying  the  theorem  and  suggests  the  next  step  to  be  taken  to  derive  an  even  more  realistic  model. 

To  apply  the  phasor  approach  of  Eq.  (2)  to  the  problem  of  modelling  the  scattered  return  from  m 
pulses,  consider  the  following  phasor  model: 

pjtsx  Ni 

E\{F,t)  =  \  <:  =  1,  .  .  .  ,  m,  (101) 

i  =  1 

where  E\  is  the  scattered  field  for  the  ifcth  of  m  pulses.  The  aspect  of  the  problem  that  should  be 
noted  is  that  in  general  each  pulse  arises  from  a  different  irtstantancous  number  of  scatterers.  To 
proceed  further  in  the  development  of  a  statistical  model  for  this  situation,  one  would  need  a  multidi¬ 
mensional  model  for  the  dynamics  of  the  number  fluctuations  from  pulse  to  pulse.  Note,  however,  that 
if  one  assumes  that  the  radar  samples  the  environment  much  more  rapidly  than  the  timescale  on  which 
the  number  of  scatterers  changes,  then  each  pulse  return  consists  of  the  same  number  of  scatterers,  and 
the  N^’s  in  Eq.  (101)  may  be  replaced  by  a  single  N.  This  situation  may  then  be  accommodated 
within  the  framework  of  the  theorem  as  follows:  Assemble  the  m  elementary  scatterer  contributions 
from  each  of  tlie  N  elementary  scatterers  into  a  random  vector  of  length  2m  (in  which  the  polar 
representation  of  the  phasor  approach  has  been  converted  to  an  inphase  and  quadrature  representation). 
The  total  return  (i.e.,  from  all  m  pulses)  then  will  consist  of  the  sum  of  these  N  random  vectors, 
which  is  precisely  the  situation  addressed  in  the  theorem.  If  this  vector  of  length  2m  is  represented  as 
a  complex  random  vector  of  length  m ,  then  the  PDF  of  the  random  vector  Y  after  application  of  the 
theorem  may  be  written 


/y(y)  = 


1 


(7i)'"I<i>I6 


/ 


1  - 

- e 

rT 


"  dF^{x), 


(102) 


where  O  represents  the  normalized  correlation  matrix  that  describes  the  pulse-to-pulse  correlation 
between  the  pulses.  (This  matrix  has  been  assumed  to  be  invertible.  For  practical  scenarios,  this 
assumption  is  justified.)  Note  that  in  the  application  of  the  theorem,  this  quantity  is  assumed  to  be 
known.  This  assumption  is  consistent  with  the  idea  of  measuring  this  quantity  for  application  to  practi¬ 
cal  problems. 

Although  the  model  described  by  Eq.  (102)  is  the  same  model  as  described  by  Conte  and  Longo, 
the  motivation  for  using  this  model  here  differs  from  the  motivation  of  Conte  and  Longo,  who  sought 


40 


NRL  REPORT  9367 


to  develop  a  model  for  the  radar  clutter  process.  They  did  seek  to  interpret  the  random  process  giving 
rise  to  such  a  multidimensional  PDF  physically  in  terms  of  a  locally  Gaussian  process  whose  power 
level  is  modulated  because  of  the  radar’s  large-scale  sampling  of  the  environment  Their  use  of  this 
model,  nonetheless,  was  somewhat  arbitrary.  Herein,  the  problem  was  approached  from  the  basic 
phenomenological  picture  of  the  scattering  process  with  minimal  assumptions  about  the  statistics  of  the 
scatterers  involved.  The  model  then  emerges  as  the  end  product  of  the  limit  theorem  presented  in  Sec¬ 
tion  2.  This  approach  complements  the  work  of  Conte  and  Longo  by  returning  to  basic  physical  prin¬ 
ciples  to  obtain  the  model;  it  is  hilly  analogous  to  the  classical  approach  of  justifying  the  multivariate 
Gaussian  model  on  the  basis  of  the  CLT.  Unlike  in  the  work  of  Conte  and  Longo,  though,  nothing  is 
stated  here  about  clutter  as  a  random  process.  The  result  here  is  limited  to  the  description  of  the  statis¬ 
tics  of  a  random  vector  of  pulse  returns.  For  this  work  to  be  extended  to  a  description  of  the  random 
process  describing  clutter  returns,  a  dynamical  model  of  the  number  fluctuations  included  in  this 
modelling  s^proach  is  needed.  Because  investigation  of  the  role  of  number  fluctuations  in  scattering 
problems  would  seem  to  be  a  new  area  of  research,  much  more  work  in  this  area  is  needed  before  an 
accurate  model  of  radar  clutter  as  a  random  process  can  be  formulated. 

6.  DETECTION  OF  TARGETS  IN  NON-GAUSSIAN  NOISE 

6.1  Structure  of  the  Optimal  Detector 

In  this  section,  the  optimal  detection  structure  for  detecting  a  signal  in  additive  noise  modelled  as 
above  is  studied.  The  resulting  structure  is  a  funcdon  of  an  optimal  estimator  of  a  random  quantity. 
This  structure  reveals  an  intimate  relationship  between  optimal  detection  and  optimal  estimation  for  this 
class  of  processes. 

Consider  the  following  hypothesis  test: 

Hq:x  =  y 

Hi  :  X  =  y  +  s  (103) 


where 


X  =  observed  m  -dimensional  complex  data  vector 


y  =  m  -dimensional  complex  noise  vector 
s  =  known  m  -dimensional  complex  signal  vector 


=  (if,  P  =  =  complex  number. 
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If  the  noise  is  modelled  as  indicated  above,  tl^n  jc  has  the  following  PDF  under  hypothesis 

i0=0,  1): 


where 


/i(jc)  = 


1 


(Tir  l<Dl  6 


<7,  =(x  -  sfO  \x  -  s),  i  =  1 
=  1=0 


(104) 


<I)  =  normalized  correlation  matrix  (assumed  to  be  invertible) 


X  =  variance  of  underlying  Gaussian  vector 


f  ^(t)  =  CDF  of  X 


I  I  =  matrix  determinant 


H  =  complex  conjugate  transpose. 


Equation  (104)  may  be  interpreted  as  the  PDF  of  a  conditional  Gaussian  random  vector  with  variance  X 
averaged  over  the  variation  in  x.  Then 


E{x\ 


(105) 


where  E  denotes  expected  value. 

The  Neyman-Pearson  likelihood  ratio  (LR)  statistic  is  given  by 


X(x)  = 


f\{x) 

/o(^) 


Define  a  fimction  as 


oo 

P„(P)  =  f  ’  dPxW 

0  XT 


(106) 


(107) 
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Then  the  LR  statistic  is  given  by 


Xix)  = 


F„i(x  -sfO  Hx  -s)) 
F^(x»0-h) 


(108) 


This  stnicture  fonns  the  starting  point  for  the  analysis  of  detection  performance  given  later  in  this 
report.  First,  however,  an  interpretation  of  the  LR  statistic  as  a  so-called  estimator-correlator  (Schwartz 
1975)  is  briefly  reviewed.  This  interpretation  is  interesting  as  it  reveals  the  LR  statistic  presented 
above  to  be  equivalent  to  a  matched  filter  compared  to  a  data-dependent  threshold. 

6.1.1  Optimal  Detector  as  an  Estimator-Correlator 

This  section  presents  a  brief  review  of  a  result  previously  given  by  Sangston  and  Gerlach  (1989). 
It  is  known  that  the  conditional  Gaussian  PDF  is  a  multivariate  member  of  the  exponential  family,  i.e., 
the  PDF  is  of  the  form 


fi(xh)  =  c(x)  h(x)  expj 


X  Tiy(x)r,(jt)} 

U  =  1 


(109)) 


where  in  this  case 


c(t)  = 


1 

(nr 


I  =  1 


T^(x)  =  sufficient  statistic  for 


T 


Reparameterize  the  conditional  Gaussian  PDF  by  setting  a  =  1/x  to  obtain 
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Define  a  function  ^  (z )  for  z  ^  0  by 


(icr  i<i>i  0 


From  Eq.  (Ill)  may  be  obtained 


ding  {2) 


(7c)'"I<I)!^(z)6 


J  dFaiOL). 


Since  the  minimum  mean-square  estimate  of  a  under  hypothesis  1  is  given  by  the  conditional  mean 
estimate  (CME)  (Poor  1988), 


Ei[Oi\x]  = - - f  dF aia) 


where  qi  is  as  defined  above,  Eq.  (112)  combined  with  Eq.  (113)  yields 


£,[abc]  =  - 


1  dgiqj) 
gi(li)  dqi 


Since 


dlngiqj) 

dqi 


Inf  j  (jc )  =  Ing  (q,)  =  |  dz  +  Ing  (0), 

n  wZ 


then,  with  the  obvious  notation  suggested  by  Eq.  (1 15),  the  structure  of  the  likelihood  ratio  is  given  by 


f  (x) 

^(^)  =  4r^  =  expU  E[a!z]dz-f  E[a\z]dz\ 
fo(x)  i  6 


where  E[a\z]  =  -dlng(z)/dz.  Equation  (116)  shows  explicitly  the  relationship  between  optimal 
detection  and  optimal  estimation  for  this  type  of  iK)ise. 
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Define  a  function  Fiq)  as 

F(q)  =  jE[a\z]dz,Q>0.  (117) 

0 

Because  £[alz]  is  nonnegative  for  positive  z  and  q  is  positive,  F(q)  is  a  monotonically  increasing 
function  of  q.  In  this  new  notation,  the  log-likelihood  ratio  may  be  rewritten  as 

Hi 

Fiqo)-F(qO^  T,  (118) 

"o 

or  equivalently, 

Fiqo)-T^  FiqO.  (119) 

"o 

Now,  since  Fiq)  is  monotonically  increasing. 

Hi 

F-\F{qQ)-T)^  qy  (120) 

Ho 

Finally,  since 

qi  =  -  2Re{s^ 0~^x)  + 

=  qQ-2Re(s^<!r^x)  + (121) 

substitution  in  Eq.  (120)  yields 

Hi 

2Re(5"^-*j:)  J  q^- F-\Fiqo) -T) -i- drh.  (122) 

Ho 

The  optimal  test  is  seen  to  be  equivalent  to  a  matched  filter  compared  to  a  data-dependent  threshold. 
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6J  Generalized  Likelihood  Ratio  Detector 


Normally  for  the  radar  detection  problem,  the  complex  signal  amplitude  of  the  desired  signal  vec¬ 
tor  is  unknown  to  the  detector.  If  unknown  parameters  are  present  in  the  LR  detector,  a  scheme 
known  as  the  generalized  likelihood  ratio  test  (GLRT)  is  commonly  used  whereby  the  unknown  param¬ 
eters  are  replaced  by  maximum  likelihood  (ML)  estimates  under  each  hypothesis  (Van  Trees  1968).  In 
this  section,  the  GLRT  for  detection  of  a  signal  with  unknown  complex  amplitude  in  additive  noise 
modelled  as  described  earlier  is  derived. 

Set  s  =  as  Qj  where  a  is  the  unknown  complex  amplitude  of  the  desired  signal  s  and  .Sq  is 
known  and  has  magnitude  one,  i.e.,  I  UqI  I  =  Sq  =  1.  The  GLRT  for  this  problem  is  thus  given  by 


X(x)  = 


F^((x  -  asQfO-\x  -  a^o)) 
F^(jc"0~‘jc) 


(123) 


where 


d  =  arg  max  F„((x  -  osq)^^  *(jc  -  osq))  (124) 

a 


is  the  maximum  likelihood  estimate  of  the  complex  amplitude  a.  As  usual,  the  maximum  of  F^  with 
respect  to  fl  can  be  found  by  setting  dF^/da  =0  and  solving  for  a.  With 
P  =  (x  -  osq^^  -  osq),  this  derivative  may  be  written  as 


^  dFm  dp 
da  d^  da’ 

From  Eq.  (107),  one  may  show 


dP 


(P)  =  -/ 


.1 


0  X 


m+l 


dF^(x). 


(125) 


(126) 


Because  F^(x)  is  nondecreasing  it  follows  easily  that 


dF. 


m 


require 


ire  4^  =  0. 


dP 


<0. 


da 


Only  one  such  solution  exists  and  is  given  by 


Thus,  the  solutions  of  interest  here 


d(x)  = 


SQd>  ^x 


(127) 


An  interesting  aspect  of  this  particular  result  is  that  it  does  not  depend  on  the  nature  of  the  CDF 
F^ix).  Substitution  of  this  result  into  the  GLRT  in  place  of  a  yields 
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XOc)  = 


-  IpP)) 
F^(jc"<D-ix) 


(128) 


where  I  I  denotes  magnitude,  and 


lpP  = 


\sq^O 

(x"<D-ix)(so"<l>-^5o)  ■ 


(129) 


It  follows  from  the  Schwartz  Inequality  that  0  ^  Ipl  ^  1. 

If  this  estimate  of  the  unknown  complex  signal  amplitude  is  substituted  into  the  matched  filter 
interpretation  of  the  LR  test,  the  result  is  easily  ^wn  to  be 

-■  <  ^0-F-\FM-T).  (130) 

•^0  ®  "^0  ir 

”0 


This  structure  is  again  seen  to  be  a  relatively  simple  stmcture  compared  to  a  data-dependent  threshold, 
where  this  same  simple  structure  is  compared  to  a  fixed  threshold  if  the  noise  is  Gaussian.  In  this 
interpretation,  the  effect  of  deviating  from  Gaussian  noise  is  to  induce  a  variability  in  the  setting  of  the 
threshold. 

63  Equivalent  GLRT 

If  O  is  nonsingular,  then  there  exists  a  nonsingular  matrix  A  such  that  the  effect  of  processing 
the  input  vector  x  by  A  is  to  whiten  the  input  rx>ise  vector,  normalize  the  input  noise  power  of  each 
element  in  the  input  noise  vector  to  one  and 


Aso  =  0, ....  0)^  =  So  (131) 

(Horn  and  Johnson  1985).  Note  that  all  of  the  desired  signal  energy  has  been  placed  into  the  first  ele¬ 
ment  by  the  matrix  transform.  Thus  if  u  =  Ax ,  then  under  Hq 

E[uu^]=I„  (132) 

where  is  the  m  X  m  identity  matrix.  Transforming  the  input  vector  by  the  nonsingular  matrix  A 
does  not  change  the  statistical  properties  of  the  GLRT  of  the  preceding  section.  Under  this  transforma¬ 
tion,  the  GLRT  becomes 


X(m)  = 


F^(m»m(1  -  Ipl^) 
Fmiu^u) 


(133) 
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where 


(u"u)(So"lo) 


lUi|2 

2  ki' 


(134) 


The  GLRT  finally  beoHnes 


Fm(2  IxJ^) 

ku) - ^ - .  (135) 

/^»(2  kP) 


Figure  S  shows  a  functional  block  diagram  of  the  equivalent  GLRT.  In  this  figure,  T  is  the  detector 
threshold. 


Fig.  5  —  Equivalent  Detector 
6.4  Performance  Assessment  of  Detector 

This  section  presents  closed-fonn  solutions  for  the  probabilities  of  detection  Pq  and  false  alarm 
Pf  associated  with  the  GLRT  presented  earlier.  The  desired  sigiuil  amplitude  is  assumed  to  fluctuate 
statistically  as  a  complex  circular  Gaussian  process  with  input  power  equal  to  (Helstrom  1968). 
As  a  result,  the  desired  signal  magnitude  is  Rayleigh  distributed  with  variance  Under  the  Hi 
hypothesis,  the  signal  after  the  A  matrix  transformation  is  contained  entirely  in  the  first  element  or 
channel  and  is  also  a  complex  circular  Gaussian  process  with  power  equal  to  =  SQ^<lr^SQ 
Note  that  at  this  point  in  the  processing  (after  the  A  matrix  transformation),  the  output  of  the  first 
channel  is  the  matched  filter  output.  The  signal-to-noise  (S/N)  ratio  associated  with  the  matched  filter, 
denoted  (5/A()^,,  can  be  expressed  as 

(S/N)^  =a,2  5o«<I)-V  (136) 
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We  now  analyze  the  statistically  equivalent  GLRT  presented  earlier.  Let 

w  =  !ir,|2 


and 


V  = 

n=2 


(137) 


(138) 


Note  that  the  random  variables  luil^,  1^2!^,  •  -  •  ,  lum  conditioned  on  X  are  statistically 

independent.  Thus,  w  and  v  when  conditioned  on  x  are  statistically  independent  Furthermore,  under 
Hq  or  //  j,  (n  =  2,  .  .  .  ,  m )  are  i.i.d.  circular  Gaussian  random  variables  with  variance  equal  to 
X.  It  is  straightforward  to  show  that  v  when  conditioned  on  X  under  Hq  and  //]  has  a  2{N  -  1)  chi- 
square  PDF  (Pa^ulis  1965)  with  =  0.5  and  is  given  by 


039) 

where  /  (  lx)  denotes  that  the  PDF  is  conditioned  on  x. 

Under  Hq,  Ui  is  a  complex  circular  Gaussian  random  variable  with  variance  equal  to  x  ,  and 
under  //j,  it  is  die  same  except  the  variance  equals  +  X.  Thus,  w  has  an  exponential  PDF  with 


/h.Ix(w  I't)  =J 


W 

1  ”7 
—e  ^ 
X 

1 


T 


0,  +T 


&,^  +  x 


Ho 

Hi. 


(140) 


The  decision  rule  in  terms  of  w  and  v  is 


Fm(y) 

Fm(y^  +V) 


Hi 


(141) 


It  was  shown  earlier  that  is  strictly  monotonically  decreasing.  As  a  result,  T  ^  1.  Since  is 
also  a  one-to-one  function,  an  equival^t  decision  rule  is 


w 


Hi 

> 


< 


-I 


Hr 


^m(v) 


-V  =  7’o. 


(142) 
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where  is  the  inverse  function  of  .  Define  P{p\v,  x)  to  be  the  probability  of  detection  when 
conditioned  on  v  and  X.  It  is  well-known  (DiFranco  and  Rubin  1980)  that  the  conditional  probability 
of  detection  of  a  random  variable  with  exponential  PDF  and  variance  ^  -I-  X  is 


/*  (Z)  Iv ,  X)  =  e 


U 

T 


e,  +T 


(143) 


From  Eqs.  (136),  (142),  and  (143),  it  follows  that 


=  j  j  /vlT(vlt)/^(x)exp 
0  0 


F^-\F„iv)  /  T)  -  V 


iS/N)^,  +  X 


dvdx. 


(144) 


where /x('t)  is  specified  and  /vix(v  It)  is  given  by  Eq.  (139). 
The  Pp  is  found  by  setting  (S/N)^,  =  0  in  Eq.  (144): 


Pp  =  i  j  /vlt(vlt)/x(t)  exp 
0  0 


F^-HF^(v)  I  T)  -  V 


dvdx. 


(145) 


6.4.1  Results 

In  this  section,  we  present  some  performance  curves  of  Pp  and  Pq  for  a  specific  Rayleigh  mix¬ 
ture  model.  The  underlying  mixing  distribution  is  chosen  to  be  the  gamma  distribution  with  PDF 
given  by 


/x(x)  = 


r(v) 


(146) 


where  r  is  the  gamma  function,  T|  determines  the  mean  of  the  distribution,  and  V  controls  the  devia¬ 
tion  from  Rayleigh  statistics.  Figure  6  shows  curves  that  describe  the  behavior  of  the  gamma  PDF  for 
different  values  of  v.  As  may  be  seen  from  the  figure,  the  statistics  of  the  Rayleigh  mixture  approach 
pure  Rayleigh  statistics  as  v  approaches  «».  Without  loss  of  generality,  Ti  is  set  equal  to  1,  which 
merely  mirnalizes  the  elemental  noise  power  to  a  specific  value.  With  this  mixing  PDF,  the  statistics 
of  the  univariate  magnitude  of  x  under  the  Hq  hypothesis,  i.e.,  the  Rayleigh  mixture  statistics,  are 
described  by  the  K  distribution 


/ixl(lj^l)  =  \x\)  (147) 

where  Ky  is  a  modified  Bessel  function  of  the  second  kind. 
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Fig.  6  —  The  gamma  PDF 


The  Pp  results  are  a  function  of  m ,  v,  and  T;th&P])  results  are  a  function  of  m ,  v,  (S /N)gpi , 
mi  Pp.  In  the  latter  case,  m,  w,  mi  Pp  are  ^cihed,  and  a  unique  threshold  T  is  determined  from 
Eq.  (14S).  This  threshold  is  then  used  in  Eq.  (144)  in  conjunction  with  m,  v,  and  (S/N)gp,  to  deter¬ 
mine  Pq  . 

For  V  =  oo,  this  corresponds  to  /t(t)  =  5(t  -  1)  where  5  is  the  Dirac  delta  function.  For  this 
case,  the  complex  input  random  variables  are  unconditionally  Gaussian,  and  it  is  known  (DiFranco  and 
Rubin  1980)  that 

Pp  = 

and 

_  r 

Pd  =  ^ 

_ 1 

_  +  1) 

Figures  7  and  8  present  curves  of  Pp  vs  In  T  for  v  =  1,  2,  10,  30,  «>,  and  m  =  2  and  5, 
respectively.  As  expected,  Pp  is  monotonically  decreasing  in  T  but  the  threshold  must  increase  to 
maintain  a  constant  Pp  for  an  increase  in  the  variance  of  X  (or  equivalently,  a  decrease  in  v.)  Also 
comparison  of  Figs.  7  and  8  show  that  for  constant  Pp  and  v,  the  threshold  decreases  for  increasing 
m. 


(148) 


(149) 
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Fig.  7  —  Pp  for  Rayleigh  mixture,  m  =^1 


Fig.  8  —  Pp  for  Rayleigh  mixture,  m  =  5 


Figures  9  and  10  present  curves  of  Pq  vs  (S//V)op,  for  v  =  1,  2,  10,  30,  oo,  Pp  =  10^,  and 
m  =  2  and  5,  respectively.  For  constant  (S/N)gpf,  Pq  decreases  as  the  variance  of  X  increases.  For 
m  =  2  and  for  a  constant  level  of  Pq  (in  the  range  0.1  to  0.9),  there  is  approximately  7  dB  difference 
in  {SlN)of^  from  v  =  1  to  v  =  oo.  For  m  =  5  and  for  a  constant  level  of  Pq  (in  the  range  0.1  to 
0.9),  there  is  approximately  4.5  dB  difference  in  (SIN)^,  from  v  =  1  to  v  =  oo. 
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Fig.  9  —  Pq  for  Rayleigh  mixture,  m  =  2,  Pp  =  10“^ 


Fig.  10  —  P[)  for  Rayleigh  mixture,  m  -  5,  Pp  =  10^ 

7.  SUMMARY  AND  CONCLUSIONS 

To  summarize,  the  problem  of  the  detection  of  targets  against  a  background  of  non-Gaussian, 
correlated  clutter  was  studied.  The  investigation  began  by  returning  to  the  underlying  phenomenologi¬ 
cal  picture  of  the  physical  scattering  process  and  deriving  the  required  statistical  model  from  that  start¬ 
ing  point  The  statistics  were  obtained  through  the  application  of  a  new  limit  theorem  from  probability 
theory  that  was  presented  and  proven  herein.  This  theorem  is  an  extension  of  the  CLT  to  the  situation 
in  which  the  number  N  of  random  vectors  being  added  together  is  itself  a  random  variable.  Since  the 
case  of  nonrandom  N  may  be  viewed  as  a  degenerate  case  of  fluctuations,  this  new  theorem  includes 
the  CLT.  This  theorem  was  then  used  to  study  tl^  problem  of  modelling  the  one-dimensional  statistics 
of  clutter  amplitude  statistics  and  to  fonnulate  a  multidimensional  model  of  clutter  statistics  that  incor¬ 
porates  correlation  between  pulse  returns.  This  multidimensional  model  formed  the  starting  point  for 
an  investigation  into  the  problem  of  radar  target  detection  against  a  background  of  non-Gaussian,  corre¬ 
lated  clutter. 
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The  woilc  presented  herein  raises  a  number  of  interesting  questions.  First,  further  justification  of 
this  type  of  analysis  to  scattering  and  propagation  problems  requires 

a.  rigorous  demonstration  that  number  fluctuations  occur, 

b.  demonstration  that  N  is  large,  and 

c.  a  methodology  for  determining  the  PMF  of  the  number 
fluctuations  from  the  i^ysics  of  the  problem. 

Limited  results  along  these  lines  are  available.  For  propagation  of  waves  through  continuous  random 
media,  Dashen  (1987)  uses  the  path  integral  formalism  to  draw  some  conclusions  about  a  and  b.  With 
respect  to  c,  Jakeman  (1980)  has  shown  that  a  birth-death-immigration  process  for  number  fluctuations 
leads  to  the  negative  binomial  distribution.  In  a  study  of  wave  propagation  through  a  turbulent  atmo¬ 
sphere,  Gochelashvili  and  Shishov  (1978)  assume  the  exponential  mixture  form  for  the  distribution  of 
intensity.  They  then  obtain  the  form  of  F,(r)  from  considerations  of  the  spectrum  of  turbulence. 
Through  the  mechanism  of  canonical  PMF’s  presented  above,  a  model  for  F,(0  may  be  used  to  obtain 
a  model  for  the  PMF  of  the  number  fluctuations.  Thus  the  work  of  Gochelashavili  and  Shishov  sug¬ 
gests  that  models  for  the  number  fluctuations  may  be  obtainable  from  knowledge  of  the  spectrum  for 
problems  of  propagation  through  random  media.  Further  work  is  needed  for  all  three  questions. 

Further  work  on  the  limit  theorem  itself  is  also  suggested.  In  particular,  one  would  like  to  know 
if  the  assumption  of  identically  distributed  random  vectors  can  be  dropped  (but  stUl  assume  indepen¬ 
dence)  or  if  the  assumption  of  independence  can  be  dropped  (but  still  assume  identically  distributed 
random  vectors.)  Both  of  the  assumptions  have  been  relaxed  in  the  case  of  the  classical  CLT,  so  the 
possibility  exists  that  they  may  also  be  relaxed  here.  These  results  would  greatly  extend  the  applicabil¬ 
ity  of  these  ideas  to  scattering  problems. 

Finally,  an  investigation  into  dynamical  models  for  number  fluctuations  is  needed  if  the  model  is 
to  be  extended  to  describe  the  radar  clutter  process,  as  opposed  to  the  description  only  of  the  statistics 
of  clutter  as  a  random  vector,  as  was  done  here. 

For  the  detection  problem,  the  interpretation  of  the  optimal  detector  in  tenns  of  the  estimator- 
correlator  structure  suggests  various  approaches  to  suboptimal  implementation  of  the  detector.  These 
problems  are  as  yet  unexplored.  Also,  the  structure  of  the  optimal  detector  with  the  proper  dynamics 
of  number  fluctuations  taken  into  account  remains  an  open  question.  Intuitively,  one  expects  the  struc¬ 
ture  to  continue  to  be  that  of  the  estimator-correlator,  although  at  this  point  this  idea  is  only  a  conjec¬ 
ture.  The  detection  structure  for  this  problem  appears  to  be  a  very  interesting  question. 

ACKNOWLEDGMENTS 

The  authors  thank  Drs.  Eric  Mokole  and  Dennis  Trizna  of  the  Naval  Research  Laboratory  for 
their  help  during  the  development  of  this  work. 

REFERENCES 

Beckmann  P.  and  A.  Spizzichino.  1987.  The  Scattering  cf  Electromagnetic  Waves  from  Rough  Surfaces 
(Artech  House,  Inc.,  Norwood  MA),  Ch.  7. 

Bernstein  S.  1929.  "Sur  les  Fonctions  Absolument  Monotones,"  Acta  Mathematica  52(1),  1. 


54 


NRL  REPORT  9367 


Boas  R.P.  1939.  "The  Stieltjes  Moment  Problem  for  Functions  of  Bounded  Variation,"  Bull.  Am. 
Math.  Soc.  45(6),  399-404. 

Bochner  S.  1937a.  "Completely  Monotone  Functions  of  the  Laplace  Operator  for  Torus  and  Sphere," 
Duke  Mathematical  Journal  3(5),  488-503. 

— .  1937b.  "Stable  Laws  of  Probability  and  Completely  Monotone  Functions,"  Duke  Mathematical 
Journal  3(9),  726-728. 

Chow  S.C.  and  H.  Teicher.  1978.  Probability  Theory:  Independence.  Interchangeability,  Martingales 
(Springer- Verlag,  New  York),  Section  9.4. 

Conte  E.  and  M.  Longo.  1987.  "Characterization  of  Radar  Gutter  as  a  Spherically  Invariant  Random 
Process,"  lEE  Proc.  Part  F  134(2),  191-197. 

Courant  R.  and  D.  Hilbert.  1953.  Methods  of  Mathematical  Physics:  Vol.  1  (John  Wiley  and  Sons, 
New  Yoik),  pp.  103-105. 

Cramer  H.  1970.  Random  Variables  and  Probability  Distributions,  3rd  ed.  (Cambridge  University 
Press,  London),  p.  103. 

Dashen  R.  1979.  "Path  Integrals  for  Waves  in  Random  Media,”  J.  Math.  Phys.  20(5),  894-920. 

— .  1984.  "Distribution  of  Intensity  in  a  Multiply  Scattering  Medium",  Opt.  Lett.  10(4),  110-112. 

Difranco  J.V.  and  W.L.  Rubin.  1980.  Radar  Detection  (Artech  House  Norwood,  MA). 

Flatte  S.M.,  S.A.  Reynolds,  and  R.  Dashen.  1987.  "Path  Integral  Treatment  of  Intensity  Behavior  for 
Rays  in  a  Sound  Channel,"  J.  Acoust.  Soc.  Am.  82(3),  967-972. 

Gochelashvili  K.S.  and  V.I.  Shishov.  1978.  "Strong  Fluctuations  of  Laser  Radiation  Intensity  in  a  Tur¬ 
bulent  Atmosphere  -  the  Distribution  Function,”  Sov.  Phys.  JETP  47(6),  1028-1030. 

Grosjean  C.C.  1989.  "An  Inverse  Laplace  Transform  -  Solution,"  SIAM  Review  31(4),  685-686. 

Helstrom  C.W.  1968.  Statistical  Theory  of  Signal  Detection.  2nd.  ed.  (Pergamon,  New  York). 

Horn  R.A.  and  C.A.  Johnson.  1985.  Matrix  Analysis  (Cambridge  Univ.  Press,  New  York). 

Hille  E.  and  J.  Tamaikin.  1933,  1934.  "On  Moment  Functions,"  Proc.  Nat.  Acad.  Sci.  19,  902-912,  20, 
140-144. 

Jakeman  E.  and  P.N.  Pusey.  1976.  "A  Model  for  Non-Rayleigh  Sea  Echo,"  IEEE  Trans.  Ant.  Prop. 
24(6),  806-814. 

— .  1978.  "Significance  of  K  Distributions  in  Scattering  Experiments,"  Phys.  Rev.  Lett.  4(9),  546-550. 
Jakeman  E.  1980.  "On  the  Statistics  of  the  K  Distribution,"  J.  Phys.  A.  13(1),  31-48. 


SANGSTON  AND  GERLACH 


Kolmogorov  AJ>I.  and  S.V.  Fomin.  1970.  Introductory  Real  Analysis  (Dover  Books,  New  Yortc). 

Olin  l.D.  1982.  "Amplitude  and  Temporal  Statistics  of  Sea  Gutter,"  lEE  Cortference  Publication  No. 
216.  Radar  82  (lEE  Press,  London). 

Pi^ulis  A.  1965.  Probability.  Random  Variables,  and  Stochastic  Processes  (McGraw-Hill,  New 
Yoik). 

Pollard  H.  1946.  "The  Representation  of  as  a  Laplace  Integral,"  Bull.  Amer.  Math.  Soc.  52,  908- 
910. 

Poor  H.V.  1988.  An  Introduction  to  Signal  Detection  and  Estimation  (Springer- Verlag,  New  Yoik). 

Pmdnikov  A.P.,  Y.A.  Brychkov,  and  0.1.  Marichev.  1986.  Integrals  and  Series:  Vol.  1  (Gordon  and 
Breach  Science  Publishers,  New  York). 

Sangston  K.J.  and  K.  Gerlach.  1989.  "Results  on  the  Etetection  of  Signals  in  Spherically  Invariant 
Noise,"  NRL  Report  9202. 

Schoenberg  l.J.  1938.  "Metric  Spaces  and  Completely  Monotonic  Functions,"  Annals  of  Mathematics 
39(4),  811-841. 

Schwartz  S.C.  1975.  "Conditional  Mean  Estimates  and  Bayesian  Hypothesis  Testing,"  IEEE  Trans.  Inf. 
Theory  21(6),  663-665. 

Skolnik  M.l.  1980.  Introduction  to  Radar  Systems.  2nd  ed.  (McGraw-Hill  Book  Company,  New  York), 
p.  496. 

Tamarkin  J.D.  1931.  "On  a  Theorem  of  S.  Bernstein- Widder,”  Trans.  Amer.  Math.  Soc.  33,  892-896. 

Trizna  D.B.  1988.  "Measurement  and  Interpretation  of  North  Atlantic  Ocean  Marine  Radar  Sea 
Scatter,"  NRL  Report  9099. 

— .  1989.  "Pacific  Marine  Radar  Sea  Scatter  Experimental  Results,”  NRL  Memorandum  Report  6506. 

Trunk  G.V.  1972.  "Radar  Properties  of  Non-Rayleigh  Sea  Clutter,"  IEEE  Trans.  Aerospace  and  Elec¬ 
tronics  Systems  8(2),  196-204. 

— 1976.  "Non-Rayleigh  Sea  Clutter:  Properties  and  Detection  of  Targets,"  NRL  Report  7986. 

Valenzuela  G.R.  and  M.  Laing.  1971.  "On  the  Statistics  of  Sea  Gutter,"  NRL  Report  7349. 

Van  Trees  H.C.  1968.  Detection,  Estimation,  and  Modulation  Theory,  Part  I  (John  Wiley  and  Sons, 
New  York). 

Ward  K.D.  1981.  "Compound  Representation  of  High  Resolution  Sea  Gutter,"  Electron.  Lett.  17(6), 
561-563. 


56 


NRL  REPORT  9367 


— 1982.  "A  Radar  Sea  Qutter  Model  and  its  Application  to  Perfonnance  Assessment,"  lEE  Confer¬ 
ence  Publication  No.  216,  Radar  82  (lEE  Press,  London). 

Watson  G.N.  1944.  A  Treatise  on  the  Theory  of  Bessel  Functions,  2nd  ed.  (Cambridge  Univ.  Press, 
New  York). 

Widder  D.V.  1931.  "Necessary  and  Sufficient  Conditions  for  the  Representation  of  a  Function  as  a 
Laplace  Integral,"  Trans.  Amer.  Math.  Soc.  33,  851-892. 

— .  1934.  "The  Inversion  of  the  Laplace  Integral  and  the  Related  Moment  Problem,"  Trans.  Amer. 
Math.  Soc.  36,  107-200. 

— .  1935.  "An  Application  of  Laguerre  Polynomials,"  Duke  Mathematical  Journal  1,  126-136. 


57 


